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Abstract 

A  polynomial  equation  for  the  eigenvalues  of  the  modes 
of  off-center  unstable  confocal  resonators  is  developed.  A 
constant  gain  for  steady  state  modes  in  a  bare  cavity  is  as¬ 
sumed.  The  field  is  built-up  from  right  and  left-traveling 
diffraction  components  for  a  number  of  round  trips  through 
the  resonator  and  geometrical  components  from  the  core  re¬ 
gion.  Using  an  asymptotic  expansion  of  the  diffraction  inte¬ 
gral,  the  boundary  conditions  are  developed.  These,  with  the 
propagation  equations  across  the  resonator,  are  used  to  re¬ 
late  the  diffraction  and  geometrical  components  to  the  diffrac¬ 
tion  amplitude  after  one  round  trip  in  the  cavity.  The  poly¬ 
nomial  equation  for  the  eigenvalues  is  developed  from  the 
first  round  trip  amplitude  function,  after  approximating  a 
slowly  varying  function  of  the  field  to  be  constant.  A  meth¬ 
od  is  proposed  for  examining  the  behavior  of  the  approximated 
function  for  the  centered  resonator  case  and  including  it  in 
mode  calculations  if  necessary. 
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ASYMPTOTIC  ANALYSIS  OF 


OFF-CENTER  UNSTABLE 
CONFOCAL  RESONATORS 

I .  Introduction 

Unstable  resonators  have  many  attractive  features  in¬ 
cluding  large  mode  volume,  substantial  transverse  mode  dis¬ 
crimination,  and  direct  output  coupling  (Ref  9:156).  Due  to 
these  characteristics,  the  determination  of  the  resonant 
modes  of  these  resonators  has  been  the  subject  of  much  re¬ 
search.  Unlike  stable  resonator  modes  (Ref  10:328-332), 
those  of  the  unstable  resonator  cannot  be  easily  described 
in  terms  of  well-known  special  functions. 

Many  methods  have  been  used  to  determine  the  modes  of 
unstable  configurations.  Geometric  approaches  (Ref  9:156) 
yield  crude  results  for  only  the  lowest-loss  mode.  These 
results  are  valid  only  for  high  Fresnel  number  resonators. 
Methods  which  iterate  to  self-consistent  field  distributions 
(Ref  2)  typically  yield  only  lowest-loss  mode  information. 

More  elaborate  techniques  have  been  developed  which  employ 
Fast-Fourier-Transform  techniques  and  include  effects  of 
gain,  diffraction,  and  index  of  refraction  variations  (Ref 
13).  This  technique,  however,  is  not  practical  for  high  Fres¬ 
nel  number  resonators  due  to  prohibitive  increases  in  the 
required  computing  time.  Other  techniques  using  numerical 


integration  also  become  costly  to  run  for  high  Fresnel  num¬ 
ber  resonators. 

An  asymptotic  theory  of  unstable  resonator  modes  using 
an  asymptotic  expansion  in  the  resonator  integral  equation 
has  been  deceloped  by  Horwitz  (Ref  3)  which  yields  the  low- 
est-loss  mode  and  higher  order  modes,  uses  little  computer 
time,  and  is  valid  for  a  wide  range  of  Fresnel  numbers.  The 
eigenvalues  of  the  modes  are  found  as  the  roots  of  a  polynom¬ 
ial  equation  in  this  method.  Horwitz  has  extended  this 
theory  to  include  misalignment  of  the  feedback  mirror  by  a 
small  tilt  (Ref  4). 

Moore  and  McCarthy  (Ref  6)  have  extended  the  asymptotic 
approach  of  Horwitz  to  determine  the  steady  state  modes  of  a 
bare  resonator  and  a  loaded  resonator  with  arbitrary  fixed 
gain  distributions.  Their  analysis  applies  to  the  positive 
branch  unstable  strip  confocal  resonator  (Ref  9:355).  For 
this  geometry  the  optical  axis  passes  through  the  center  of 
both  mirrors.  This  centered  confocal  resonator  geometry  is 
chosen  because  the  beam  is  not  focused  within  the  cavity,  al¬ 
so  the  output  beam  should  be  well  collimated  since  in  the  geo¬ 
metric  limit  it  would  be  a  plane  wave. 

In  comparison  to  the  centered  resonator  geometry,  reson¬ 
ators  with  an  off-center  feedback  mirror  appear  to  offer  ad¬ 
vantage  in  beam  quality.  Better  mode-loss  separation  and 
more  far  field  energy  in  the  first  Airy  square  (for  rectan¬ 
gular  mirror  cavities)  may  be  achieved  (Refs  7:2159  and  14: 
1828).  To  examine  these  characteristics,  a  method  of 


2 


determining  the  modes  for  off-center  unstable  confocal  res¬ 
onators  is  needed. 

The  primary  objective  of  this  work  is  to  extend  the 
analysis  of  Moore  and  McCarthy  to  obtain  u  polynomial  eq¬ 
uation  for  the  eigenvalues  of  the  off-center  geometry.  A 
secondary  objective  is  to  examine  the  approximations  made  by 
Moore  and  McCarthy  in  their  analysis  of  the  centered  geometry 
and  to  report  the  results  of  this  work. 

This  report  consists  of  six  chapters  und  five  appendices. 
Following  this  introduction.  Chapter  II  presents  the  basic 
formalism,  assumptions,  und  general  approach  used  in  this  an¬ 
alysis.  Chapter  III  develops  the  equations  for  the  boundary 
conditions  on  the  mirrors  while  the  amplitude  functions  for 
the  field  in  the  resonator  are  developed  in  Chapter  IV.  The 
polynomial  equation  for  the  eigenvalues  of  the  off-center  un¬ 
stable  resonator  modes  is  developed  in  Chapter  V.  An  approx¬ 
imation  made  by  Moore  and  McCarthy  to  obtain  the  polynomial 
equation  for  the  centered  case  is  also  examined  in  this  chap¬ 
ter.  The  last  chapter  presents  specific  conclusions  and  rec¬ 
ommendations.  The  appendices  include  mathematical  derivations 
and  detailed  discussions  supporting  the  work  presented  in  the 
main  body  of  this  report. 
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II .  Problem  Formulation 


The  purpose  of  this  chapter  is  to  set  the  framework  for 
use  in  the  following  development.  Four  major  areas  are  dis¬ 
cussed.  The  first  area  is  the  basic  formalism  for  the  prob¬ 
lem.  This  includes  a  discussion  of  the  geometry,  assumptions, 
and  general  approach  used  in  the  development.  The  second 
section  presents  the  basic  form  assumed  for  the  field  in  the 
cavity.  Next,  gain  in  the  cavity  is  considered  in  light  of 
assumptions  made  for  steady  state  modes  to  be  supported. 
Finally,  the  assumed  form  of  the  field  amplitude  expressions 
is  discussed. 

Basic  Formalism 

Geometry.  The  geometry  of  the  off-center  unstable  strip 
confocal  resonator  considered  in  this  report  is  illustrated 
in  Fig.  (1).  This  cavity  consists  of  a  convex  cylindrical 
feedback  mirror  at  the  right  end  and  a  large  concave  cylin¬ 
drical  mirror  at  the  left  end.  The  feedback  mirror,  M^,  is 
located  distance  d  to  the  left  of  the  common  focal  point, 
f.p.  The  large  mirror  is  displaced  a  distance  D  to  the 
left  of  the  feedback  mirror  along  the  optical  axis.  The  op¬ 
tical  axis  is  taken  as  the  line  joining  the  centers  of  curv¬ 
ature  of  the  mirrors,  C ^  and  C2  . 

The  upper  edge  of  the  feedback  mirror  is  located  dis¬ 
tance  b  from  the  optical  axis  and  the  lower  edge  distance 


Fig  1  OFF-CENTER  UNSTABLE  CONFOCAL  GEOMETRY 


-a.  For  this  development  it  is  assumed  that  |b|>ja|.  The 
transverse  dimensions  of  the  large  mirror  are  assumed  to  be 
sufficiently  large  so  that  all  of  the  incident  radiation  is 
reflected  without  introduction  of  diffraction  effects.  Prac¬ 
tically,  this  "sufficiently  large"  transverse  dimension  can 
be  considered  to  be  at  least  b(M+l)  (Ref  6:228)  where  M 
is  the  resonator  magnification  developed  in  Appendix  B  and 
given  by 


M 


D+d 

~r 


(i) 


Two  effective  Fresnel  numbers  can  be  defined  for  this 
cavity  due  to  the  asymmetry  of  the  feedback  mirror  edges  to 
the  optical  axis.  They  are 

Feb  "  TO  (2) 

Fea  "  JSa  <3> 

These  are  also  reffered  to  as  equivalent  Fresnel  numbers  ' 
(Refs  3:1529  and  12:39). 

The  asymmetry  caused  by  the  feedback  mirror  not  being 
centered  on  the  optical  axis  can  be  introduced  by  starting 
with  the  centered  geometry  and  rotating  the  feedback  mirror 
about  its  center  of  curvature.  The  difference  between  the 
asymmetry  being  introduced  by  a  rotation  rather  than  a  tilt¬ 
ing  of  the  feedback  mirror  is  discussed  in  Appendix  A.  It 
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is  clear  that  when  the  feedback  mirror  edge  distances  from 
the  optical  axis  are  equal,  |b|*|a|  ,  the  expressions  de¬ 

veloped  in  this  analysis  should  reduce  to  those  of  the  cen¬ 
tered  case. 

As  in  Moore  and  McCarthy 's  analysis  (Ref  6:228),  the 
field  in  the  cavity  can  be  modeled  by  a  right-traveling  mod¬ 
ulated  plane  wave  and  a  left-traveling  modulated  cylindrical 
wave  originating  at  the  focal  point.  To  describe  these  fields 
conveniently,  two  coordinate  systems  are  used,  see  Fig.  (1). 

A  cartesian  coordinate  system  (z,x)  with  its  origin  at  the  in¬ 
tersection  of  the  optical  axis  with  the  large  mirror,  has  pos¬ 
itive  z  defined  to  the  right  of  the  mirror  and  positive  x 
above  the  axis.  The  cylindrical  coordinate  system  (p,0)  has 
its  origin  at  the  focal  point.  Positive  0  is  defined  as  a 
clockwise  rotation  in  the  (z,x)  plane  about  the  focal  point 
with  0*0  along  the  optical  axis  to  the  left  of  the  focal 
point. 

When  the  paraxial  approximation  is  valid  these  coordin¬ 
ate  systems  are  related  by 

p  ■  Md-z 

0  -  i 

P 

Then  the  paraxial  coordinates  (p,x)  can  be  used  to  describe 
locations  relative  to  the  focal  point. 

Assumptions.  To  make  the  development  more  tractable, 


(4) 

(5) 
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without  sacrificing  accuracy,  five  major  assumptions  are  made: 

1.  The  field  is  scaler. 

2.  The  mirrors  are  perfectly  conducting. 

3.  The  laser  is  operating  in  steady  state. 

4.  The  main  contributions  to  the  modes  are  from  paraxial 
regions  in  the  resonator. 

5.  The  optical  axis  does  not  pass  closer  than  several 
Fresnel  zones  from  the  feedback  mirror  edge. 

The  first  assumption  greatly  simplifies  the  analysis  of 
the  problem.  The  scalar  equations  neglect  the  vector  nature 
of  the  field  and  so  polarization  effects  are  not  modeled. 

This  allows  the  use  of  the  strip  resonator  geometry.  In  a 
bare  cavity  analysis  this  is  not  a  serious  drawback.  If  a 
specific  gain  medium  was  being  modeled  then  polarization  be¬ 
comes  a  factor. 

The  assumption  of  perfectly  conducting  mirrors  implies 
that  all  of  the  incident  light  is  reflected.  Thus  losses  due 
to  absorption  by  the  mirrors  are  not  modeled. 

Oscillation  in  steady  state  at  a  single  frequency  neglects 
time-dependent  phenomena.  It  also  implies  a  constant  gain  in 
the  resonator  that  exactly  cancels  losses  due  to  output  coup¬ 
ling.  This  analysis  differs  from  that  of  a  truly  bare  cavity 
(no  gain  medium)  which  would  exhibit  decaying  modes. 

The  assumption  of  paraxial  modes  implies  that  the  major 
contribution  to  the  mode  comes  from  a  paraxial  region  about 
the  geometric  source  point. 

8 

,  ....  . . 


It  is  assumed  that  the  optical  axis  does  not  pass  closer 
than  several  Fresnel  zones  from  the  edge  of  the  feedback  mir¬ 
ror  because  of  the  asymptotic  expansion  used  in  this  analysis. 
This  expansion  is  carried  out  to  terms  of  the  order  of  the  in¬ 
verse  Fresnel  number  for  the  edge  nearest  the  optical  axis. 

If  the  optical  axis  passes  closer  than  several  Fresnel  zones 
this  expansion  becomes  increasingly  inaccurate. 

Approach.  The  general  approach  to  obtaining  the  poly¬ 
nomial  equation  for  the  mode  eigenvalues  consists  of  four 
parts.  The  first  part  assumes  a  form  for  the  field  in  the 
cavity  built  up  from  a  geometric  contribution  and  diffraction 
components  due  to  the  feedback  mirror  edges.  The  boundary 
conditions  relating  these  field  compenents  are  developed  in 
the  second  part  of  the  analysis.  In  this  section  the  devel¬ 
opment  of  the  boundary  conditions  on  the  feedback  mirror  in¬ 
volves  the  use  of  an  asymptotic  approximation  for  a  diffrac¬ 
tion  integral  evaluated  over  the  limits  of  the  mirror.  The 
third  part  of  the  analysis  develops  expressions  relating  the 
geometrical  components  and  diffraction  components  of  the  am¬ 
plitude  function  to  the  diffraction  components  after  one 
round  trip  in  the  resonator.  The  diff’ action  amplitude  com¬ 
ponents  are  then  expressed  in  terms  of  an  auxiliary  slowly 
varying  function.  In  the  fourth  part  an  approximation  of  the 
auxiliary  function,  similar  to  Moore  and  McCarthy's  (Ref6:233), 
is  made  to  obtain  a  polynomial  for  the  eigenvalues.  This  sec¬ 
tion  then  returns  to  the  centered  case  and  proposes  a  matrix 
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equation  that  could  be  used  to  examine  the  auxiliary  function 
approximation  and  determine  its  validity. 

Field  Considerations 

Field  in  Resonator.  The  field  inside  the  cavity  is  as¬ 
sumed  to  consist  of  a  right-traveling  wave  and  left-traveling 
wave  similar  in  form  to  those  of  geometrical  theories  (Ref  9: 
157).  The  right-traveling  wave  has  the  form  of  a  modulated 
plane  wave  where  its  amplitude  function  includes  rapidly  os¬ 
cillating  transverse  phase  terms  due  to  diffraction  effects. 
Similarly,  the  left-traveling  wave  assumes  the  form  of  a  mod¬ 
ulated  cylindrical  wave  originating  at  the  focal  point  of  the 
cavity.  Its  amplitude  function  also  includes  diffraction  ef¬ 
fects.  The  total  field  in  the  cavity  is  then  given  by 

E  ■  |f(z,x)exp{ik(z+Md)}-p~ig(p,9)exp{ikp}Je~iwt  (6) 

The  mirrors  of  the  cavity  impose  boundary  conditions  up¬ 
on  the  field.  At  the  mirrors  the  phases  of  the  right  and  left 
traveling  waves  must  match.  Also,  the  total  field  must  vanish 
on  the  two  mirrors.  This  yields  the  boundary  conditions 

f (0,M6d)  -  (Md)~*g(Md,9)  (7) 

d*f(D,9d)e2ikD;-  |<9<§ 
g(d,9)  -  d  a  (8) 

0  jotherwise 
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Gain  Considerations 


From  the  assumptions  used  in  this  analysis,  the  field 
must  obey  the  scalar  Helmholtz  equation 


V2E 


92  1 

3t2  c*" 


E 


2 yia  9 
c  9t 


(9) 


where 

c  ■  velocity  of  light 
U  =  permeability  of  gain  medium 
o  =  conductivity  of  gain  medium 
Vio  =  G  the  gain  coefficient  of  the  medium 
Since  the  field  consists  of  right  and  left-traveling  waves, 
these  traveling  waves  individually  satisfy  the  Helmholtz 
equation.  Substituting  the  right-traveling  wave  from  Eq.(6) 
into  Eq.(9)  yields 

2ik^f(z,x)-Gf(z,x)^  +  1^-  f(z,x)  *=  0  (10) 

when  second  derivatives  of  the  z-coordinate  are  neglected  due 
to  the  slow  variation  of  the  field  as  a  function  of  z. 

Similarly,  when  the  left-traveling  wave  of  Eq.(6)  is  sub¬ 
stituted  into  Eq.(9) 

2ik^g(p,e)-Gg(p,9)^  +p"2|jrg(P,9)  “  0  (ID 

where  second  derivatives  with  respect  to  p  are  neglected  due 
to  the  slow  variation  of  the  field  as  a  function  of  p  .  Al¬ 
so  terms  of  order  ip”^g(p,6)  are  neglected  when  compared  to 
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2kp  ^  ~  g(p,6)  .  In  Eq.(lO)  and  Eq.(ll)  the  G  term  is  as- 

<7  P 

sumed  to  be  a  constant. 

If  diffraction  effects  are  ignored  the  second  derivative 
terms  of  Eq.  (10)  and  Eq.(ll)  can  be  neglected  (Ref  6:229) 
yielding  the  geometrical  rate  equations 


|^f(z)  «  Gf (z)  (12) 

!^g(p)  ®  Gg( p )  (13) 

The  solutions  of  Eq.(12)  and  Eq.(13)  are 

f(z)  **  f(0)exp{Gz)  (14) 

g(p)  *  g(0)exp{Gp>  (15) 


The  right  and  left-traveling  waves  in  the  cavity  are  de¬ 
pendent  upon  each  other.  This  dependence  can  be  expressed  in 
terms  of  the  gain  in  the  cavity.  The  Wronskian  of  Eq.(12)  and 
Eq.(13)  can  be  found  by  multiplying  Eq.(12)  by  g(p)  and  Eq. 
(13)  by  f(z)  and  then  subtracting  one  equation  from  the  oth¬ 
er.  In  the  paraxial  approximation,  the  p  dependence  is  chang¬ 
ed  to  z  using  Eq.(4).  This  yields 

fcH  g(Md-z)  -  |ig(Md-z)jf(z)  (16) 

Integrating  yields 

f(z)g(Md-z)  -  C  (17) 
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where  C  is  a  constant  of  integration. 

Evaluating  this  expression  at  each  mirror  yields 

f (O)g(Md)  -  f(D)g(d)  (18) 

Upon  substituting  for  f(0)  from  Eq.(7)  and  for  f(D)  from 
Eq.(8)  it  is  found  that 

e2lkDM-ig2(Md)  -  g2(d)  (19) 

From  Eq . ( 15 ) 

g2 (Md)  =  g2(d)e2GD  (20) 

A  constant  gain  coefficient  Q  which  satisfies  the  steady 
state  requirement  for  this  analysis  can  be  defined  as 

Q  »  GD 

Then  substituting  Eq.(20)  into  Eq.(19)  yields 

2ikD  2Q 
U  «  M  ae  e  x 

where,  for  this  geometrical  case 

U  -  1  (23) 

The  y  terms  are  the  eigenvalue  for  the  resonator  modes. 
When  diffraction  effects  are  included,  the  u  terms  are 
found  to  be  the  roots  of  a  polynomial  equation  where  each 
root  is  the  eigenvalue  of  a  different  resonant  mode.  Thus 


(21) 


(22) 
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the  constant  gain  required  to  sustain  each  mode  in  steady 
state  is  different. 


Diffraction  Effects 

The  amplitude  functions  of  the  right  and  left-traveling 
waves  are  expanded  in  terms  of  a  geometrical  amplitude  func¬ 
tion  and  edge  diffraction  amplitude  functions  as  shown  by 
Horwitz  (Ref  3:1530)  and  Moore  and  McCarthy  (Ref  6:230).  The 
diffraction  terms  consist  of  slowly  varying  amplitude  functions 
and  phase  terms  that  have  rapid  transverse  oscillations.  Well 
within  the  geometrical  shadow  boundaries  the  diffraction  terms 
can  be  viewed  as  cylindrical  wavelets  emanating  from  the  vir¬ 
tual  images  of  the  feedback  mirror  edges  (Ref  6:230).  From 
Appendix  B,  these  virtual  images  are  located  on  two  half¬ 
parabolas.  For  the  (d,b)  edge  of  the  feedback  mirror,  the 
images  are  located  at 

(pn»bn)  =  (dM2n,bMn)  ;n=0, ±1, ±2 , . . . , ±N  (24) 

and  for  the  (d,-a)  edge 

(pn*-an)  “  (dM2n,-aMn);n=0,±l,±2,...,±N  (25) 

For  n>0  these  images  lie  behind  the  large  mirror  and 
their  cylindrical  wavelets  contribute  to  the  right-traveling 
wave  amplitude  function  which  can  be  written  as 
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f(z,x)  = 


N 


,(pn-p)  3rn(p,x)expj 


U=1 


ik  (x-bn  )‘ 

5-  r^r 


+(Pn-p)-isn(p,x)exp 


ik  (x+an  )2 I 
2  (Pn-P)  j 


+f(p,x) 


(26) 


where 


rn(p,x)  =  diffraction  amplitude  functions  from  left 
images  of  edge  (d,b) 

sn(p,x)  =  diffraction  amplitude  functions  from  left 
images  of  edge  (d,-a) 

A 

f(p,x)  =  geometrical  amplitude  function 
The  n=0  wavelets  originate  at  the  feedback  mirror  edges 
and  the  n<0  wavelets  originate  at  the  virtual  images  to  the 
right  of  the  feedback  mirror.  These  images  approach  the  focal 
point  as  n-<— °°  .  These  cylindrical  wavelets  contribute  to  the 

left-traveling  amplitude  function  which  takes  the  form 


p"*g(p.e) 


-N+l 

Y  (p-pn)"*un(p,x)exp 

ik 

2~ 

n=0,^l,-2 

1  1 

L 

pn  P 

(27) 


+p  2g(p,X) 
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where 


un(P,x)  ■  diffraction  amplitude  function  from  right 
images  of  edge  (d,b) 

vn(p,x)  =  diffraction  amplitude  function  from  right 
images  of  edge  (d,-a) 

A 

g(p,x)  ■  geometrical  amplitude  function 
The  phases  appearing  in  Eq.(26)  and  Eq.(27)  are  the  Fres¬ 
nel  approximations  to  the  phases  of  cylindrical  waves  emanat¬ 
ing  from  points  (Pn»bn>  or  (Pn,-an)  '  When  lnl  ls  sma11 
the  phase  terms  are  rapidly  oscillating  but  when  |n|  becomes 

large  the  phases  associated  with  edge  (d,b)  terms  approach 
the  constant  phase 

exp{^H  =  exp{2TriFeb}  (28) 

and  the  phases  associated  with  edge  (d,-a)  terms  approach 
the  constant  phase 

exp,rtcr}  =  exp{2lTiFea}  (29) 

It  has  been  found  by  Horwitz  (Ref  3:1533)  that  a  practical 
value  for  large  N  is  given  by 

MN  =  250Feff  (30) 

This  is  practical  in  the  sense  that  a  larger  N  changes  the 
eigenvalues  at  most  in  the  third  decimal  place.  In  the 
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r 
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off-center  case  let  F  „  -F  ,  since  F  ,  '•F  .  This  will  ns- 

eff  eb  eb  on 

sure  a  large  enough  number  of  terms  in  the  expansion. 

Finally,  substituting  Eq.(2G)  and  Eq.(27)  into  Eq.(9) 
yields  the  rate  equations 


3rn 


3s 


Gr 


JL  - 


Gs. 


3? 

T57 
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(31) 


Gun 


“  Gv 


n 


A 


C.g 


(32) 


3 

whore  is  the  directional  derivative  alon£  the  direction 

u  s 

of  the  diffraction  ray  propagation  off  the  optical  axis.  This 
is  consistent  with  previous  gain  considerations. 

In  the  special  case  where  |a|  ■  |b|  these  preceding  ex¬ 
pressions  reduce  to  those  of  Moore  and  McCarthy  (Ref  6)  for 
the  centered  unstable  con  focal  resonator. 

In  summary,  this  section  has  developed  the  expansions  for 


the  right  and  left-traveling  wave  amplitude  functions  in  terms 
of  geometric  and  edge  diffraction  amplitude  functions. 


III.  Boundary  Conditions  on  the  Mirrors 


This  chapter  details  the  development  of  the  boundary  con¬ 
ditions,  including  the  diffraction  terms,  on  the  two  mirrors 
of  the  cavity.  These  conditions  specify  the  relationships  be¬ 
tween  the  right  and  left-traveling  waves.  The  development  is 
arranged  into  two  sections.  The  boundary  conditions  on  the 
large  mirror  are  developed  first.  The  second  section  devel¬ 
ops  the  boundary  conditions  on  the  feedback  mirror  using  an 
asymptotic  expansion  of  a  diffraction  integral. 

Boundary  Conditions  on  the  Large  Mirror 

Since  the  large  mirror  is  assumed  to  reflect  all  of  the 
incident  field,  the  boundary  conditions  on  it  are  easily  ob¬ 
tained.  From  Eq.(7)  the  boundary  conditions  on  the  large  mir¬ 
ror  are  found  be  equating  the  right  hand  sides  of  Eqs.(26)  and 


(27)  for  p=Md  and  0“x/Md  .  Using  Eq.(24)  it  can  be  shown 
that  the  resulting  equation  has  like  exponentials  on  each  side 
of  the  equality 


The  coefficients  of  these  like  exponentials  must  also  be  equal. 
Equating  these  coefficients  yields  the  boundary  conditions  for 
the  diffraction  amplitudes: 

(Pn-Md)-irn(Md,x)  =  (Md-P1_n)"4u1_n(Md,x)  (35) 

;n«-l,2, .  . .  ,N 

(pn-Md)"isn(Md,x)  =  (Md-P1_n)_4v1_n(Md,x)  (36) 

;n*l,2, . . . ,N 

Equating  the  remaining  terms  yields  the  boundary  condition  for 
the  geometrical  amplitudes: 

f ( Md , x )  =  (Md)~*g(Md,x)  (37) 

Boundary  Conditions  on  the  Feedback  Mirror 

The  boundary  conditions  on  the  feedback  mirror  are  more 
difficult  to  develop.  While  the  total  field  must  vanish  on 
the  mirror,  the  left-traveling  component  must  be  zero  off  the 
mirror.  A  diffraction  integral  evaluated  over  the  feedback 
mirror  will  account  for  these  requirements.  It  is  assumed  that 
the  gain  coefficient  is  zero  very  near  the  feedback  mirror,  so 
this  term  is  not  included  in  the  integration.  Then,  the  Huy- 
gen-Fresnel  diffraction  integral  will  give  an  expression  for 
g(p,0)  in  terms  of  g(d,6): 


10 


g(p,e) 


g(d.e') 


(38) 


where  the  Fresnel  degree  of  approximation  has  been  used  in 
phase  terms.  Making  the  change  of  variables,  —  , 

and  replacing  g(d,0')  with  the  form  from  Eq.(8)  yields 


The  expression  for  f(D,an)  given  by  Eq.(26)  is  substituted 
into  this  equation.  This  integral  can  be  approximated  by  an 
asymptotic  expansion  about  the  stationary  phase  point  as  de¬ 
scribed  in  Appendix  C.  The  results  of  the  specific  integra¬ 
tion  in  Eq.(39)  are  developed  in  Appendix  D  where  the  approx¬ 
imation  neglects  terms  of  higher  order  than  F  “l  .  The 
expression  for  g(p,0)  near  the  feedback  mirror  is  then 
found  to  be 
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where  Ij,I2,  and  Ij  are  developed  in  Appendix  D.  Eq.(D— i: 
of  this  appendix  presents  Eq.(40)  in  detail. 

Recalling  Eq.(27)  and  letting  the  summation  index  n  = 


yields 

'N-l 

g  (P,9)  =  p*(p-p-m)"* 

_m=0, 1 , 2 


The  expressions  for  g(p,0)  in  Eq.(41)  and  Eq.(40)  can  then 
be  equated.  It  can  be  shown  that  as  p  approaches  d  ther 
are  like  exponentials  on  both  sides  of  the  equality 


n  ■  1,2. . . , N-l 
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and 


Therefore  the  coefficients  of  these  like  exponentials  can  be 
equated  as  p  approaches  d  .  These  are  the  boundary  con¬ 
ditions  for  the  diffraction  amplitude  functions  on  the  feed¬ 
back  mirror: 


2ikD,  N  .  , 

(pn"d)  3rn(d,x)  =  (d-p_n 


)~*u_n(d,x) 


(44) 


n  =  1,2, .. . ,N-1 


2ikD .  ,.-i  .  ,  . 

e  (Pn“d)  sn(d»x> 


(d-p_n)"*v-n(d»x) 


(45) 


n  =  1,2, .. . ,N-1 


The  geometrical  term  from  Eq.(41)  can  then  be  equated  to  the 
geometrical  term  and  the  Nth  diffraction  terms.  Here  N  is 
assumed  large  as  given  by  Eq.(30)  and  so  the  phases  are  con¬ 
stant  as  given  by  Eq.(28)  and  Eq.(29).  The  boundary  condi¬ 
tion  for  the  geometrical  terms  as  p  approaches  d  is  then 


d'^d.x)  -  e21kD[(pN-d)-1rN(d,x)exp{i||!}  (4g) 

+(PN-d)“isN(d,x)exp|i||-}+f(d,x)] 
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Then  the  m«0  terms  remain  from  Eq.(41)  which  equal  the  re¬ 
maining  terms  from  Eq.(40).  This  is  taken  as  the  expression 
for  the  diffraction  amplitude  functions  on  their  first  trip 
across  the  cavity  and  is  used  in  the  next  chapter. 

In  summary,  this  chapter  has  developed  the  relationship 
between  the  various  amplitude  functions  at  the  mirror  bound¬ 
aries.  It  remains  to  relate  these  amplitude  functions  to  a 
common  expression,  the  diffraction  amplitudes  after  one  round 
trip  in  the  cavity. 
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IV.  Amplitude  Functions  in  the  Resonator 

This  chapter  develops  the  relationships  among  the  geo¬ 
metric  and  the  diffraction  amplitude  functions.  The  devel¬ 
opment  of  these  relationships  consists  of  three  sections. 

The  first  section  relates  the  diffraction  amplitude  functions 
on  the  nth  round  trip  of  the  cavity  (equivalent  to  emanating 
from  the  nth  virtual  image  of  the  feedback  mirror  edges)  to 
the  first  round  trip  value.  The  first  round  trip  amplitudes 
are  then  defined  in  terms  of  auxiliary  functions,  q  (x)  and 
qfa(x)  .  The  second  section  relates  the  geometrical  amplitude 
functions  to  these  auxiliary  functions.  In  the  last  section 
the  expressions  for  these  auxiliary  functions  are  developed 
in  terms  of  the  eigenvalues  for  the  resonator  modes. 

It  is  important  to  determine  these  relationships  be¬ 
tween  the  various  amplitude  functions  because  the  values  of 
these  amplitudes  can  be  found  once  the  value  for  the  first 
round  trip  diffraction  amplitude  functions  are  known.  With 
values  for  all  the  amplitude  functions  the  total  amplitude 
functions  for  the  right  and  left-traveling  wave  can  be  deter¬ 
mined. 

Diffraction  Amplitude  Functions 

Initial  Amplitude  Functions.  After  the  boundary  condi¬ 
tions  on  the  feedback  mirror  are  extracted  from  the  right 
hand  side  of  Eq.(41)  equated  with  the  right  hand  side  of 
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Eq.(40),  the  remaining  terms  yield  expressions  for  the  diffrac¬ 
tion  amplitude  functions  on  their  first  trip  across  the  cavity 
from  the  feedback  mirror.  From  these  remaining  terms  it  can 
be  shown  that  the  following  exponentials  are  equal  by  using 
Eq.(24)  and  Eq.(25): 


/ 


exp 


ik 

2~ 


This  leads  to  the  expressions  for  the  diffraction  amplitude 
functions  on  their  first  trip  across  the  cavity: 


Uo(P,x) 
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v0(p,x) 


"  ♦  (ot) 


2ikD 


N 


■<(Vd)'4rn<d’-a) 


n«l 


x  exp 


»ik  (a+bn Ta+bn  k  e(a+0d) 

\^~  r^-Sr (  +  d(p-d) 


-1 


(50) 


x“>»-d>~Vd'-‘>oxP{ir  ^3^! 


-1 


+  f ( d, -a ) 


P  ( a+Od) 


d( p-d) 


-1 


The  +  at  the  beginning  of  the  righthand  side  of  Eqs.(49)  and 
(50)  explicitly  denotes  the  choice  of  roots  from  the  square 
root.  In  Eq.(40)  the  positive  root  of  the  square  root  was  im¬ 
plicitly  chosen  so  that  the  real  part  was  positive  (Ref6:231). 
This  leads  to  the  choice  of  the  negative  root  here. 

In  Eq.(49)  it  can  be  shown  that 


p(b-0d) 

d(p-d) 


-1 


and 


where 


(51) 


(52) 


(53) 
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is  the  angle  with  respect  to  the  optical  axis  of  the  ray  join¬ 
ing  the  mirror  edge  (d,b)  with  a  point  (p,x)  in  the  cavity, 

*bn  "  +  ^  "  |(l+M“n)“1  (54) 

is  the  angle  for  the  ray  joining  image  point  (p  ^ , b  )  with 
edge  (d,b), 


(58) 


(59) 
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is  the  angle  of  the  ray  joining  edge  (d,-a)  with  a  point 
(p,x)  in  the  cavity, 


is  the  angle  of  the  ray  joining  points  (P_n»b_n)  with  edge 
(d,-a), 


(61) 


is  the  angle  of  the  line  joining  points  (p_n<-a  n)  with 
edge  (d,-a)  .  These  angles  are  illustrated  in  Fig. (3).  A1 

so  it  can  be  shown  that 
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Fig  3  RAY  ANGLES  FOR  v0(p,x)  AMPLITUDE  FUNCTION 


p(a+9d) 
d(p-d) 

It  can  also  be  shown  that  using  the  angle  expressions,  a 
and  8  ,  allows  the  left-hand  sides  of  Eq.(49)  and  Eq.(50) 

to  be  written  as 


U0(p,x)  =  u0(p, (P— d)6+b) 

(63) 

v  0 ( P , x )  =  v0(p,(p-d)a-a) 

(64) 

Recall  that  Eq.(49)  and  Eq.(50)  are  valid  only  near  the  feed¬ 
back  mirror  because  the  gain  coefficient  was  assumed  to  be 
zero  when  these  equations  were  developed.  To  propagate  ac¬ 
ross  the  cavity,  maintaining  steady  state,  the  diffraction 


-1 


-i-l 


3  * 


(62) 
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amplitudes  of  Eq.(49)  and  Eq.(50)  must  be  multiplied  by  the  ex 
ponenti&l  of  the  gain  coefficient,  Q  ,  from  Eq.(21).  Then, 
using  Eq.(51)  through  Eq.(64),  the  diffracted  amplitude  func¬ 
tions  on  the  first  trip  across  the  resonator  can  be  express¬ 
ed: 

u0(p,(p-d)B+b)  =  -  eQe2ikD 

N 

x  L(0„-<i)'4rn(d,b) 
n=l 

*  Hr  [ebn-*1] ‘65> 


and 


-Kpn-d)-*sn(d,b)exp{^  ■■} 

x[®an"p]”1 

+f(d,b)[j-B  1 

v  (p.(p-d)a-a)  «  -{pkfl)*  eQe21kD 
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x  2 (pn_d)”*rn(d’~a) 
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+  (Pn-dr*s  (d,.-H)oxpJ^ 

.L  .  h1  -  .  r.  i-il 


k\n+'f  +f(d,-a) 
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Similarly,  for  the  nth  ray  paths  to  the  left,  intersecting 
the  feedback  mirror  at  any  point  (d,x)  ,  yields 


un(p,x+ 


x-b 

d-p 


a(p-d))= 

n 


n  = 


un(d,x)e^ 

-1,-2 . -N+l 


(69) 


vn(p,x+f=fI1'(p"d))=  vn(d,x)eQ 


(70) 


n  -  -1,-2, .. . , -N+l 


Amplitude  Functions  from  nth  Transit  of  Cavity.  The  ex¬ 
pressions  relating  rn  to  Tj  and  to  s^  will  now  be 

developed.  The  diffracted  amplitude  function  rn(d,x)  is  the 
amplitude  of  the  diffracted  wave  caused  by  the  feedback  mir¬ 
ror  edge  (d,b)  after  it  has  made  n  round  trips  through  the 
cavity.  This  function  can  be  related  to  the  rn_j(d,x') 
plitude  function  by  propagating  it  back  through  the  cavity 
one  round  trip.  Here  (d,x')  is  the  point  a  ray  must  start 
from  on  the  feedback  mirror  to  propagate  to  point  (d,x) 

Starting  with  rn(d,x)  and  using  Eq.(67)  to  propagate 
the  amplitude  function  on  the  feedback  mirror  to  the  large  mir¬ 
ror  yields 


rn(d,x)  ■  rn(Md,n)eQ 

t 

n  -  1,2, ... ,N 


(71) 
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where,  using  Eq.(l)  and  Eq.(24) 


n  ■  x 


(72) 


At  the  large  mirror  the  boundary  condition  of  Eq.(35)  gives 


e«rn(Md,n)  -  e« 


(73) 


n  ■  1,2, ... ,N 


From  this  Eq.(69)  is  used  to  propagate  u^  n(Md,n)  back  to 


the  feedback  mirror: 


,Q  /  M2P-M  \ 

\M-M2(1-n)J  Ul-n(Md,n) 


„2Q/Mjn-M  \  ,  . 

6 


(74) 


where 


2,3, ... ,N 


(75) 


The  boundary  condition  at  the  feedback  mirror  Eq.(44)  then 


yields 


2QfM*n-M 


,2(l-n. 


ul-n(d’x' > 
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(76) 


^e^Vid.x') 


/l-Ma(l~n)  Mln-M 


n  =  2,3, ... ,N 

Using  Eq.(72)  and  Eq.(75),  the  x'  value  can  be 
terms  of  x  .  Then  the  rn  amplitude  function 
ed  to  the  r  function  using  Eq.(71),  Eq.(73), 
Eq. (76) : 


r  (d,x)  =  M4e2Qe21kD 


/.  MIn-M 

x  r  -  (  d,x— — — 

n_1  \  M2n-1 


+  bMtV(M-M~1) 
M2n-1 


) 


n  =  2,3, ... ,N 

From  this  equation  r  ,  can  be  related  to  r 

n-l  n- 

then  also  be  related  back  to  rn-  Continuing  in 
leads  to  an  expression  relating  rQ  to  r^: 

n-l 

r  (d,x)  -  m-5"  2<n-DQe21kD(n-l) 


x  r^d.^Cx))  ;  n=l,2,...N 


where 


kbn<x)  yir^ 


expressed  in 
can  be  relat- 
Eq.(74)  and 


(77) 


which  can 
this  manner 

(78) 


(79) 
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By  a  similar  procedure  the  diffracted  amplitude  function 
sR  can  be  related  to  s^  which  yields 


Sn(d*x)  =M  T1  e2(n-l)Qe2ikD(n-l) 


(80) 


X  sx(d,kan(x))  ;  n=l ,  2 ,  . . .  N 


where 


Hjf  2  _1  \  U2n-V 


k  (x)  -  X  _a  M__=M 

an  M  n-l  M2n-1 


(81) 


First  Round  Trip  Amplitude  Function.  From  Eq.(78)  and 
Eq.(80)  the  nth  round  trip  diffracted  amplitude  functions  are 
related  to  that  of  the  first  round  trip.  Expressions  are  now 
developed  for  the  first  round  trip  amplitudes. 

From  Eq.(65)  the  diffracted  amplitude  u0  is  known  at 
any  point  (Md,x')  on  the  large  mirror  by  letting  p**Md  and 
substituting  x’  into  Eq.(53)  for  angle  6  .  This  gives 

the  amplitude  function 


u0(Md,x')  =  u  0 ( Md , b+ ( Md-d ) fl ) 


(82) 


where 


„  _  x'-b 
6  d(M-l) 


(83) 


The  boundary  condition  at  the  large  mirror  Eq.(35)  is  then 
used  to  relate  u0  to  r^ 


i 


uo (Md,x’ ) 


(84) 


J 

rj (Md,x'  ) 


Then  Eq.(67)  is  used  to  propagate  rj(Md,x')  across  the  cav¬ 
ity  to  r»(d,x): 

e^ri(Md,x')  =  ri(d,x)  (85) 

Therefore,  combining  Eq.(84)  and  Eq.(85)  yields 

ri(d,x)  =  M^e^u0  (Md,  x' )  (86) 

where 

x'  =  U+b)^  (87) 

Also,  from  Eq.(83) 


xM~ 1 -bM~2 
d(l-M"2) 


(88) 


Finally,  using  Eq.(86),  Eq.(88),  Eq.(54),  and  Eq.(55)  in  Eq.(65) 
yields  the  needed  diffraction  amplitude  expression,  Tj(d,x)  , 


r, (d,x) 


2Q  2ikD 
M3e  ^e 


(89) 


N 

(1-Mn)2 

^  ^d~^(M"n-l)~^exp 

ikb2 

2d 

M2n-1 

n=l 

X  rfd.b)!*  ]■» 

[d  1-M“2n  d( 1-M”2 )  J 
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+<r*(M2n-l)-iexp 


ika2 

2d 


(t^ 

M2n-1 


(89) 


Xsn(d,b) 


^+M-n 
a  a 


xM_1-bM~2 


3  1-M-2  n  d( 1-M-2 ) 


-l 


+f(d 


UN[b  xM-1-bM“2l 

■b)l* "  ^P7J 


- 1 


Similarly,  the  amplitude  function  Sj (d,x)  can  be  dev¬ 
eloped: 


sx  (d,x)  =  -M^e2%2ikD 


N 


y^d~*(M2n-l)~*exp*  ikb 


n=l 


2d  M2D-1 


Xrn(d,-a) 


b  t+M_n 

d 


xM-1+aM-2 


1-M 


-2n 


d(l-M“2) 


-l 


+d~*(M2n-l)“*exp 


ika2 

2d 


( 1-Mn ) 2 
M2n-1 


(90) 


Xsn(d,-a) 


a 

d 


1-M~n 

1-M"2n 


+ 


xM~  ^aM"2 
d(l-M“2) 


+?(d, 


+ 


xM_1+aM~2  ] 

d(l-M-2)  J 


38 


Up  to  this  point  the  diffraction  amplitude  functions, 
uo(p,0)  and  vo(p,0)  ,  have  been  propagated  through  one 

round  trip  in  the  cavity.  Also,  all  of  the  succeeding  round 
trip  amplitude  functions  have  been  expressed  in  terms  of 
these  first  round  trip  amplitudes  (Eqs.78  and  80).  Now  at¬ 
tention  is  turned  to  finding  the  value  of  these  first  round 
trip  amplitude  functions  by  using  an  auxiliary  function  q(x) 
(Ref6:232). 

Auxiliary  Function  q(x)  From  Eq.(89)  a  new  function, 

qb(x)  ,  can  be  introduced  which  accounts  for  the  slowly  vary¬ 
ing  x-dependence  of  the  diffracted  amplitude  function  r,(d,x). 
Let 

r, (d,x)  *=  M*e2Qe2ikDqb(x)  [b-  g]“‘  (91) 


Similarly,  from  Eq.(90)  a  function,  q  (x)  ,  associated 

A 

with  diffracted  amplitude  Sj(d,x)  is  introduced: 

s i ( d , x )  *=  M4e2Qe2ikDqa(x)  [a+  g]  *  (92) 


Now  the  rn(d,b)  amplitude  function  of  Eq.(89)  can  be  ex¬ 
pressed  in  terms  of  qb(x)  .  From  Eq.(78)  at  x«b 

n-1 


rn(d,b)  -  M~^  e2(n-l)Qe2ikD(n-l) 


(93) 


X  r ,  (  d ,  k.  (  b  ) ) 


where,  from  Eq.(79) 
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bn' 


From  Eq.(91) 


bn' 


.  b  M~ 1 +M 1 ~n 

1+M~n 

-  Mie2Qe2ikD 

X  -1  (h 

b  1"M"2 

M  i+m-»  ) 

i+M~n 

(94) 


(95) 


->  -i 


Then,  substituting  Eq.(95)  into  Eq.(93) 
rn(d,b)  -  e2nQe2nikD 


X  q, 


/  r, 

\  1+M~n  /  L 


1-M 


-2 


1+M 


-n 


-l 


(96) 


Similarly,  using  Eq.(80),  Eq.(81),  and  Eq.(92) 

n 


sn(d,b)  =  M“  e2nQe2nikD 
X  q 


a\ 

[a  (l«;  M -°)  (l-H-*)l 

L  l-M-in  J 


+  a  M“n| 


(97) 


Now  an  expression  for  qfe(x)  is  needed.  Using  Eq.(91) 
and  substituting  Eq.(96)  and  Eq.(97)  into  Eq.(89)  gives  the 


following  expression  for  qfe(x) 
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approximating  Mn/(M2n-l)  ”  by  M-n^2  ,  the  q.  (x)  ex¬ 


pression  of  Eq.(98)  becomes 


q  (x)  =  V 

\2*kb2)  1“ 


id_\J  Vll  e2nQe2nikD 


X  exp 


2d  M2n 


q,  (b 

-1  I  \  1-M 


)d-M~ 

zn 


X( 1-M-2 n )  ( l-M-n )  -  bM-n(l-M-2)(l-M-n)^  -l 

d-M"  n)  (b-  g) 


~T  <^2nQ(12nikD<avn  ika2  (a  +Mn) 


+M  ^  e  ^e  exp 


(99) 


2d  M2n-1 


X  qa  b 


[1_n(l+  E  M~n)-M“‘(b  +M"n 


X  (1-H 


/a  +  M~n\  +  bM~n(  1-M~2  )  (b  +M~n) 
_  ’  ( 1-M-2 n )  (b  -  g) 


+f(d,b)dib(l-M“2) 


Finally,  a  similar  procedure  leads  to  an  expression  for 


qa<*>  : 


•  -m  is 


I  Vm't  e2nQe2nikD 


(I  +Mn) 


„  _  ikb2  +Mr 

X  eXP  T 
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X  q 


b 


M*~n(l+  a  M~n)  -M~*  (  a  +  M~n) 


1-M 


-2  n 


) 


X  ( 1-M-2  n ) 


(! 


a'!~n(  1-M~2  ) (a  +M~n) 
(1-M"2n)  (a+  §) 


-n 


...  —  2  nQ  2nikD 

+M  2  e  exp 


ika2  ( 1-Mn )2 


2d  M2n-1 


(100) 


a  \  l-M-2  n  / 


X  ( 1-M-2  n ) 


(l-M"") 


aM~n(l-M~2 )(1-M~n)2 

(1-M_2n)  ( a+  g) 


+f(d,-a)d*a(l-M-2 ) 


The  analysis  to  this  point  has  expressed  the  rR  and 
sn  diffraction  amplitude  functions  of  Eq.(26),  at  the  feed¬ 
back/output  mirror  plane,  in  terms  of  the  first  round  trip 
diffraction  amplitude  functions  rt  and  s(  .  These  first 
round  trip  amplitudes  have  then  been  expressed  in  terms  of 
the  auxiliary  functions  q^(x)  and  qa(x)  •  Using  the  re¬ 
lationships  between  the  various  diffraction  amplitudes  on 
succeeding  trip  through  the  resonator  to  the  first  round 
trip  amplitudes,  the  q(x)  functions  can  be  expressed  in 
terms  of  a  series  expansion  in  terms  of  q(x)  at  distinct 
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locations  on  the  feedback  mirror  and  the  geometrical  com¬ 
ponent,  f  ,  of  the  field  amplitude  from  Eq.(26).  Attention 
is  now  turned  to  expressing  this  geometrical  amplitude  in 
terms  of  the  q(x)  functions. 

Geometrical  Amplitude  Functions 

Expressions  for  the  geometrical  amplitude  functions 
f(d,b)  and  f(d,-a)  from  Eq.(89)  and  Eq.(90)  are  needed  to 
evaluate  the  field  amplitude.  From  the  boundary  condition 
Eq.(37)  the  geometrical  terms  can  be  propagated  to  the  feed¬ 
back  mirror  along  the  geometrical  rays. 

From  Eq. (37) 

f ( Md , Mx )  =  (Md)~*g(Md,Mx)  (101) 

The  f  wave  propagates  along  a  ray  parallel  to  the  optical 
axis.  The  f(Md,Mx)  amplitude  is  multiplied  by  the  steady 
state  gain  in  propagating  across  the  cavity  to  the  feedback 
mirror.  Then 

f ( d , Mx )  =  eQf(Md,Mx)  (102) 

The  g(Md,Mx)  amplitude  of  Eq.(lOl)  must  be  propagated  back¬ 
wards  to  its  origin  on  the  feedback  mirror.  This  ray  origin¬ 
ates  at  the  focal  point  and  intersects  the  feedback  mirror  at 
(d,x)  .  The  steady  state  gain  must  be  included  in  this  prop¬ 

agation  expression  also: 

g(Md,Mx)  =  eQg(d,x) 


(103) 


Then  from  Eq.(102)  and  En.(103) 


?(d,Mx)  *  M~^e2<^d_^g(d,x) 


(104) 


Recalling  the  boundary  condition  Eq.(46)  at  the  feedback  mir¬ 
ror,  this  expression  can  be  written  as 


f ( d , Mx ) 


2Q  2ikD 

M  J  e  e 


X 


(pN-d)"*rN(d,x)exp 


ikb2 

TcT 


(105) 


+  (pN~d)~*SN(d,Jc)e?cp  ^3" 


+f (d,x) 


From  this  expression  for  rN  and  sN  in  terms  of  qb(x)  and 
q  (x)  are  needed.  From  Eq.(78) 

u 

rNvd,X>  -  e2(N-l )Qe21kD(N-l ) 

( 106 ) 

X  Md-kbn(x» 

where,  for  N  large 

kb„(x)  -  §  (107) 
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From  Eq. (91) 


ri(d,kbN(x))  *=  M4e2Qe2ikD 


X<*b(si)  bd-M-1) 


Then,  combining  Eq.(108)  and  Eq.(106) 


rN(d,x)  =  M2  e2NVNikD 


A  % 


(m) 


b( 1-M-2  ) 


(108) 


(109) 


Similarly,  the  expression  for  the  s  term  is  found  to 

n 


sN(d,x)  =  M2  e  we 


7  «2NQ  2NikD 


>i"n) 


a(l-M~2 ) 


(HO) 


Then,  approximating  MN/(M2N-1)  4  by  M-N/2  ,  using 

U  from  Eq . (22 ) ,  and  the  right  hand  sides  of  Eq.(109)  and 
Eq.(llO)  in  Eq.(105)  yields 

f(d,Mx)  =  yld"4uN 


X  qb(fi)  exp  ^§j-||l>(l-M_’) 

x  qa(“S)  exp  [a(1~M~2> 


+f(d,x) 
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The  geometrical  expression,  f(p,x)  ,  is  the  amplitude 
of  a  plane  wave  moving  toward  the  right.  In  the  purely  geo¬ 
metrical  case  where  the  left  mirror  is  "sufficiently  large" 
to  reflect  all  the  incident  wave,  this  amplitude  is  constant 
across  the  wave-front.  Then 


J(d,Mx)  =  f(d,x) 


(112) 


From  this,  Eq.(lll)  yields 

N+l 


f (d,x)  =  d_i  H 


1-y 


-i 


b( 1-M~* ) 


_ 

-i ' 

\  1  ika2 ) 

,  -  W-2  V 

)exp  1  2d  1 

a(l-M  ) 

1 

4 

(113) 


Noting  that  this  expression  has  no  x-dependence,  as  expected, 
it  then  follows  that 


f(d,b)  =  f (d,-a)  =  f ( d , x ) 


(114) 


Auxiliary  Functions  in  Polynomial  Equations 

The  qb(x)  and  q&(x)  auxiliary  functions  of  the  rj(d,x) 
and  Sj(d,x)  amplitude  functions  can  now  be  written  in  terms 
of  a  polynomial  expression  for  the  eigenvalues.  Using  Eq.(113), 
Eq.(114),  Eq.(2),  Eq.(3),  and  Eq.(22)  in  Eq.(99)  gives 
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r 


Similarly, 
qa(x)  - 


for  Qa(x) 


'  2*  fc) 


N 


n=l 


exp 


2iriF 


(I  +Mn) 


eb  MIn-l 


n 


X  qb(  b 


0 


M~ 1  (l+  b  M~n)  -M*~n  (E  +M~n) 


1-M 


-2n 


X  (1-M"2n) 


/b  +M-n\  +  aM~n(l-M~2)(a  +M~n) 

'  ^  /i  _  .  X\ 


-1 


(l-M  (a  +  £) 


+  exp  2  it  iF 

r  An  n  n 


ea  M2n-1 


n 


(116) 


( 1-M~n ) ( M~ 1 +M 1 
1-M‘2n 


X  ( l-M-2  n  ) 


(1-M"n) 


aM~n( 1-M~2 ) ( 1-M~n ) 2 
(1-M_2n)  (a  +  £) 


+ 


+ 


qb(n)  exp|2”1Febl 

' 

exp  |  2xiFea]  1  ■ 
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In  summary,  the  rn  and  sn  diffraction  amplitude  func¬ 
tions  of  Eq. (26)  were  related  to  the  amplitudes  rt  and  Sj 
These  first  round  trip  diffraction  amplitude  functions  were  re¬ 
lated  to  the  auxiliary  functions  qfe(x)  and  qa(x)  •  The  geo¬ 
metrical  amplitude  function,  f  ,  from  Eq.(26)  was  also  re¬ 
lated  to  the  q.  ( x )  and  q  (x)  functions.  Equations  (115)and 
(116)  are  the  resulting  expressions  for  the  auxiliary  functions. 
The  next  chapter,  using  several  assumptions,  develops  the  eigen¬ 
value  polynomial  from  these  equations.  Each  eigenvalue  can 
then  be  used  to  solve  for  the  rn  *  sn  ’  and  *  of  E(l*(26) 
for  each  mode. 
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Eigenvalue  Equations  and  Limitations 


The  purpose  of  this  chapter  is  the  development  of  ex¬ 
pressions  from  which  the  eigenvalues  of  the  off-center  un¬ 
stable  resonator  can  be  determined  and  their  limitations. 

The  first  section  develops  a  polynomial  equation  using  an  ap¬ 
proximation  on  the  q(x)  functions.  The  eigenvalues  are  the 
roots  of  this  equation.  The  second  section  discussed  the 
validity  of  the  approximations  made  in  the  first  section. 

It  then  proposes  a  matrix  equation  which  can  be  used  to  de¬ 
termine  the  eigenvalues  and  the  behavior  of  the  q(x)  function. 
The  final  section  discusses  the  limitations  on  the  application 
of  the  eigenvalue  expressions. 

Polynomial  Expression  Assuming  Constant  q(x) 

Development  of  Polynomial.  To  obtain  a  tractable  poly¬ 
nomial  equation  for  the  eigenvalues,  jj  ,  of  the  centered 
resonator,  Moore  and  McCarthy  (Ref  6:233)  make  two  approxi¬ 
mations.  First,  the  function  q(x)  is  assumed  to  be  a  con¬ 
stant.  Secondly,  with  the  exception  of  exponential  terms, 

M-n  is  neglected  when  compared  to  one.  To  obtain  the  equiv¬ 
alent  expressions  for  the  off-center  case  requires  the  approx- 
imat ions 


qb(x>  -  % 


(117) 
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(118 


(X)  -  q 


-n  I  a  b  I 

M  <<minimum  of  |l,F,or  a| 


(119! 


Here,  qb  and  q&  are  constants.  In  the  case  of  the  center 
ed  resonator  q&=qb  •  Equation  (119)  does  not  apply  for  ar 
guments  of  exponential  terms.  Using  these  approximations  in 
Eq.(115),  and  dividing  through  by  qb  yields 


*■  KO 


12  u”  exP  2"iFeb 

n=l  M  -1 


+  |  exp|2iriF 

qb  a  ea  M2n-1 


PN+1 

r —  exp 


2iriF 


qa  b 

+  q^  a  exp  2,iF€ 


(120 


Similarly,  after  dividing  through  by  q  rather  than  q. 

U  D 

Eq.(116)  gives 


+  exp 


(121) 


Then  Eq.(120)  and  Eq.(121)  have  two  unknown  parameters,  y 
and  (^a/Q^,)  »  which  can  be  determined.  After  multiplying 

through  Eq.(121)  by  (bq^aq^)  and  subtracting  from  Eq.(120), 
the  resulting  equation  can  be  solved  for  (q  /qK)  ,  result- 
ing  in  the  following  expression 
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N 


2>n 

n=l 


exp 


2niF  ,  .<i=“Y 
eb  MJn-l 


-  exp 


2n  iF 


(E  +M°)  ’ 


eb  Min-1 


(122) 


*  |  *  (1+  {2hb)  '  (^) 


N 


n=l 


/  b  \ 2 

exp 

2tt  iF  la  ±»Q 

ea  Man_1 

-  exp 


2*  if  U=i£ii 

ea  M!n-1 


The  eigenvalue  polynomial  in  terms  of  u  can  then  be 
obtained  by  substituting  this  expression  into  Eq.(120)  for 

qa^qb  *  Bringing  •  a11  terms  on  the  left  to  the  right  hand 
side  of  the  equation,  multiplying  through  by  (1-y)  ,  and 

collecting  like  summations  yields  the  eigenvalue  polynomial: 


°  -  EZ 

n=  i  k= l 


N 


(123) 


n*=  i 


+yN+1D+y-l 
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where 
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For  B 
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2tt1F 


(1-Mn)2 
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exp 


2iriF 
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-exp 
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-exp 


2TriF 


(I  +Mn) 
ea  M2"-! 


exp 


2u  iF  . 

eb  M2k-1 


-exp  j  2iriF 


(§  +Mk) 
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exp 
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exp 
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eb  2  n 
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(?  +Mn) 


ea  „2n 


-1 
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(125) 


The  eigenvalues  for  the  modes  of  the  off-center  unstable 
confocal  resonator  are  the  roots  of  the  polynomial  in  Eq.(123). 
In  the  case  of  the  centered  resonator  where  a=b  ,  these  ex¬ 
pressions  reduce  to  Moore  and  McCarthy's  (Ref  6:233). 

Field  on  the  Feedback  Mirror.  The  right-traveling  field 
on  the  feedback  mirror  can  now  be  determined  to  within  a  con¬ 
stant  multiplicative  factor,  q&  1  or  qb  1  .  From  Eq.(91) 

and  Eq. (92) ,  using  Eq.(22),  rt(d,x)  and  Sj(d,x)  can  be  ex¬ 
pressed  as 

r , ( d, x)  =  Myqb[b-  g]"1  (128) 

sx(d,x)  =  Myqa[a+  (129) 


Substituting  these  into  Eq.(78)  and  Eq.(79)  yields 

-  i 


rQ(d,x)  =  Mnpnqt 


b- 


kbn(x) 

M 


(130) 


- 1 


sn(d,x)  =  Mnunqa 


a+ 


kan<x> 


(131) 


Using  Eq. ( 130)  and  Eq.(131)  in  Eq.(26)  and  also  approximating 
M2n/(M2n-l)  ^  as  unity,  the  right-traveling  field  amplitude, 
to  within  the  multiplicative  constant  qfe  ,  is  given  by 

N 

f(d,x)qb_l  =  yVV 
n=l 
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Validity  of  Approximations  and  Proposed  Solution 

For  the  centered  resonator  case  the  feedback  mirror 
edges  are  an  equal  distance  from  the  optical  axis.  Then 
a=b  and  Eq.(115)  becomes 


X 


M~ 1 +M 1 
1+M“n 


X 


aM~n( 1-M~2  ) 

( 1+M~n )  ( a-  g). 


(133) 
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This  expression  reduces  to  Moore  and  McCarthy’s  (Hof  6:232) 
expression  only  if  q(x)  ■  q(-x)  for  the  case  of  symmetric 
modes  in  the  centered  geometry. 

The  approximation  that  M-n  can  be  neglected  when  com¬ 
pared  to  one  is  valid  for  large  M,  and  n  greater  than  one 
or  two.  However,  for  cases  where  the  magnification  approach¬ 
es  unity  it  can  easily  be  seen  that  neglecting  the  M~n  terms 
(for  small  n)  when  comparing  this  term  to  one  is  a  poor  approx¬ 
imation. 

In  Moore  and  McCarthy's  analysis  for  the  centered  reson¬ 
ator  geometry  the  auxiliary  function,  q(x)  ,  is  assumed  to 
bo  a  constant.  This  allows  the  q-terms  to  be  divided  out  of 
Kq.(133)  to  obtain  the  eigenvalue  polynomial.  It  is  not  clear 
whether  or  not  this  is  a  valid  approximation.  Moore  and 
McCarthy  offer  no  Justification  for  it.  A  method  of  solving 
for  the  eigenvalues  and  also  examining  the  behavior  of  the 
q(x)  function  without  neglecting  M  n  terms  can  be  develop¬ 
ed  however . 
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Equation  (133)  expresses  q(x)  in  terms  of  a  summation 
over  q(xk)  values  where  xfc  are  specific  locations  on  the 
output  mirror.  Each  q(x^)  can  be  expressed  in  the  form  of 
Eq.(133).  For  any  upper  bound,  N  ,  of  the  summation  this 
leads  to  a  system  of  2N+1  equations  of  the  form  given  in 
Eq.(133).  This  system  of  equations  can  be  written  as  a  non¬ 
linear  matrix  eigenvalue  problem.  If  this  can  be  solved  for 
the  true  eigenvalues  (ie.  no  constant  q(x)  assumption  made) 
then  the  resulting  eigenvectors  can  be  used  to  determine  the 
behavior  of  the  q(x)  function  for  various  modes  in  the  cen¬ 
tered  resonator  geometry. 

To  obtain  the  matrix  equation  form  the  following  symbol¬ 
ism  is  introduced: 


c  ‘ '  k  fa) 


(134) 


Ei0  -  exp  2„iFe 


(135) 


E#  -  exp  |  2iriFe  | 
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where  n»l,2,...,N  and  k=0,l,...N. 

Replacing  q(x)  by  q(x^)  in  the  left  hand  side  of  Eq. 
(133)  and  using  the  above  symbolism  leads  to  an  equation  for 
each  q(x^)  °f  the  form 

N 

0  *  T,  CE+nA+n<V‘><’‘+n)“n 

n=l 


+CE-nA-n(xk)q(x-n)Ml1  <140) 

N+l 

+  iqr  2CE»q(x.> 


-q<V 

The  2N+1  equations  of  this  form  (from  the  2N+1  values  of 
k  )  can  then  be  written  in  the  homogeneous  matrix  equation 
form 

A(u)*£“0  (141) 

where  A (u)  is  a  square  matrix  of  dimension  (2N+1)  and  £ 
a  vector  of  (2N+1)  elements.  The  matrix  elements  of  A(p) 
consist  of  the  coefficients  of  q(xR)  in  Eq.(140).  Appen¬ 
dix  E  develops  representative  elements  in  this  matrix.  These 
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coefficients  contain  only  one  unknown  parameter,  y  ,  the 
eigenvalue.  This  is  then  a  non-linear  matrix  eigenvalue 
problem.  Problems  of  this  type  are  very  difficult  to  solve 
in  general.  If  this  one  can  be  solved  it  will  yield  the 
true  eigenvalues,  y  ,  and  the  elements  of  the  eigenvector 
£  ,  for  each  eigenvalue  can  be  used  in  the  right-hand  side 

of  Eq.(133)  to  determine  the  behavior  of  the  function  q(x). 

Limitations 

There  are  limitations  to  the  range  of  applicability  of 
these  expressions  which  are  due  to  the  geometry  of  the  reson¬ 
ator  and  approximations  used  in  this  development. 

The  polynomial  expression,  used  to  determine  the  eigen¬ 
values  for  the  modes,  is  valid  only  for  the  positive  branch 
unstable  confocal  geometry  with  an  off-center  (or  centered  in 
the  case  where  a=b  )  feedback  mirror.  Since  a  magnification 
of  unity  cannot  be  achieved  by  this  configuration,  the  expres¬ 
sions  developed  become  increasingly  inaccurate  as  M  approach¬ 
es  one.  At  M  equal  to  one  the  expressions  blow  up  as  illus¬ 
trated  by  the  geometrical  contribution  to  the  field  in  Eq.(132). 

The  asymptotic  expansion  used  in  this  development  has  a 
singularity  at  the  geometrical  shadow  boundary.  The  eigen¬ 
value  expressions  are  correct  but  the  field  can  only  be  deter¬ 
mined  well  within  the  shadow  boundaries. 

As  an  effective  Fresnel  number  for  the  resonator  approaches 
unity,  this  analysis  becomes  increasingly  inaccurate.  This 
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stems  from  the  asymptotic  expansion  where  terms  of  higher 

order  than  F  “l  are  neglected.  As  the  effective  Fresnel 
ea 

number  approaches  unity,  more  terms  must  be  included  in  the 
series  expansion  for  the  diffraction  integral.  For  the  off- 
center  geometry  this  requires  that  the  optical  axis  not  pass 
very  close  to  an  edge  of  the  feedback  mirror  in  this  analysis. 

Another  limitation  to  this  analysis  is  due  to  the  Fresnel 
approximation  used  in  the  phase  terms  of  the  diffraction  am¬ 
plitude  functions.  When  these  amplitudes  are  used  in  the  dif¬ 
fraction  integral  the  resulting  field  amplitude  becomes  valid 
a  small  distance  from  the  feedback  mirror.  How  close  this 
distance  is  to  the  mirror  is  limited  by  the  Fresnel  approx¬ 
imation.  This  implies  that  the  left-traveling  field  at  a 
point  close  to  the  feedback  mirror  consists  of  the  contribu¬ 
tions  at  the  point  where  the  geometric  ray  intersects  the 
feedback  mirror  and  from  points  in  some  small  region  about 
that  point  that  satisfy  the  Fresnel  approximation.  Contri¬ 
butions  from  other  points  outside  this  region  are  thus  as¬ 
sumed  to  be  negligible. 


VI .  Conclusions  and  Recommendations 


The  primary  objective  of  this  work  has  been  accomplish¬ 
ed.  An  eigenvalue  polynomial  for  the  modes  of  an  off-center 
unstable  confocal  resonator  has  been  developed.  In  the 
special  case  of  the  centered  feedback  mirror,  the  off-center 
eigenvalue  polynomial  should  reduce  to  the  form  developed  by 
Moore  and  McCarthy (Ref  6:233).  This  has  been  verified,  val¬ 
idating  the  off-center  eigenvalue  polynomial  in  this  limit. 

In  fulfillment  of  the  secondary  objective,  a  non-linear 
matrix  eigenvalue  problem  has  been  developed  to  examine  the 
behavior  of  the  q(x)  function  for  the  centered  case.  If 
solutions  to  this  problem  are  feasible,  the  approximation  of 
q(x)  as  a  constant  can  be  examined. 

It  is  recommended  that  methods  of  solving  the  q(x)  non¬ 
linear  matrix  eigenvalue  problem  be  investigated.  If  solving 
the  problem  proves  feasible,  the  behavior  of  the  q(x)  func¬ 
tion  should  be  investigated.  If  q(x)  varies  significantly 
from  a  constant  value,  its  effects  on  the  modes  of  the  reson¬ 
ator  should  bf  determined. 

If  q(x)  can  be  evaluated  for  the  centered  case,  it  is 
recommended  that  this  technique  be  extended  to  examine  the  be¬ 
havior  of  the  functions  qb(x)  and  qft(x)  for  the  off-center 
case.  Their  combined  effect  in  the  off-center  case  may  make 
modes  for  this  resonator  geometry  more  sensitive  to  deviations 
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of  the  q(x)  functions  from  constant  values. 


I 


With  the  eigenvalue  polynomial  for  the  off-center  geom¬ 
etry,  the  computer  program  developed  by  Moore  and  McCarthy 
should  be  modified  to  calculate  the  eigenvalues  and  the  modes 
for  the  off-center  geometry.  A  parametric  study  could  then 
be  undertaken  of  the  eigenvalues  and  modes  of  various  reson¬ 
ator  configurations. 

This  asymptotic  analysis  should  be  extended  to  allow 
the  determination  of  the  field  beyond  the  geometrical  shadow 
boundaries.  This  could  be  accomplished  using  either  the 
method  employed  by  Moore  and  McCarthy  (Ref  6:234)  or  the  uni¬ 
form  asymptotic  expansion  of  Horwitz  (Ref  4:172).  With  the 
field  across  the  entire  output  plane  the  mode  intensity  dis¬ 
tribution  could  be  determined.  Then  Fast  Fourier  Transform 
techniques  could  be  used  to  calculate  far  field  intensity  dis¬ 
tributions. 
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Appendix  A 

Phase  Effects  of  Tilt 
Versus  Rotation  of  the  Feedback  Mirror 

The  purpose  of  this  appendix  is  to  illustrate  the  differ¬ 
ences  that  can  occur  in  the  distance  a  ray  from  the  feedback 
mirror  travels  in  an  off-center  resonator  due  to  tilting  the 
feedback  mirror  rather  than  rotating  it  about  its  center  of 
curvature.  The  difference  in  distance  traveled  is  related 
to  phase  differences. 

Rotated  Feedback  Mirror 

If  the  feedback  mirror  is  centered  on  the  optical  axis 
and  then  rotated  about  its  center  of  curvature,  an  off-cen¬ 
ter  resonator  of  the  same  geometry  as  used  throughout  this 
report  is  generated.  This  rotated  off-center  unstable  con- 
focal  resonator  is  illustrated  in  Fig. (5). 

In  this  case  the  (p,0)  coordinate  system  has  its  origin 
at  the  center  of  curvature  of  the  feedback  mirror  rather  than 
the  focal  point.  The  (z,x)  coordinate  system  has  its  ori¬ 
gin  at  the  intersection  of  the  optical  axis  and  the  large  mir¬ 
ror.  The  mirrors  are  separated  by  a  distance  D  .  The  ra¬ 
dius  of  curvature  of  the  feedback  mirror  is  2d.  Also,  the 
distance  separating  the  centers  of  curvature  of  the  two  mirrors 
is  S  .  For  the  confocal  case  S  equals  D. 
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Fig.  5  ROTATED  MIRROR  OFF-CENTER  RESONATOR 

For  ease  of  computation,  the  case  where  D  equals  d  , 
i.e.  resonators  with  a  magnification  of  two  are  considered. 

For  a  point  (2d,0)  on  the  mirror,  the  distance,  L  , 
from  it  to  the  point  where  the  optical  axis  intersects  the 
large  mirror  is  given  by 

L  =  (z2+x2 )*  ( A-l ) 

where 


z  =  D+2d(l-cos0) 


x  ■  2d  sin0 
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In  the  paraxial  range  of  0  ,  sin0=0  and  cos0=l  .  Then 

Eq.(A-l)  becomes 

L  =  (D2+4d202)*  ( A-2 ) 

In  comparison  to  the  centered  case  the  distance  is  iden¬ 
tical  for  a  similar  point  (2df0)  .  So,  no  phase  changes 

must  be  accounted  for  in  extending  the  analysis  to  the  off- 
centered  confocal  resonator  when  the  off-center  condition  is 
introduced  by  a  rotation  of  the  feedback  mirror  about  its 
center  of  curvature. 

Tilted  Feedback  Mirror 

If  the  centered  feedback  mirror  is  now  tilted  by  an 
angle  6  about  its  center,  the  optical  axis  shifts  and  an 
off-center  resonator  is  generated.  This  geometry  is  depict¬ 
ed  in  Fig. (6). 

This  resonator  is  no  longer  confocal  since  S'  ,  the 
new  separation  of  the  centers  of  curvature,  and  the  mirror 
separation,  D'  ,  are  no  longer  equal.  From  the  law  of 
cosines,  the  distance  S'  is  given  by 

S'  =  [4d2  +  (2d+S)2-4d(2d+S)cos.$|*  (A-3) 

Recalling,  that  S=d 

S'  *  d[l2(l-cosS)+l]*  ( A-4 ) 
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Fig.  6  TILTED  MIRROR  OFF-CENTER  RESONATOR 

Then  the  change  in  the  center  of  curvature  separation  (and 
also  in  mirror  separation)  is  given  by 

S'-S  =AS=[l2(l-cos6)+lj *d-d  (A-5) 

If  it  is  assumed  that  6  is  small  then  (1-cosS)  is  also 
small.  Using  the  binomial  approximation  of  the  bracketed 
term  in  Eq.(A-5) 

AS  =  6d( l-cos6 )  ( A-6 ) 

Using  the  series  expansion  of  cos  6 

AS  =  3d62  ( A-7) 
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Then  the  mirror  separation,  D'  ,  is  given  by 

D'  =  D-AS  (A-8) 

Now  the  distance,  L'  ,  to  feedback  mirror  point  (2d,0) 
from  the  intersection  of  the  optical  axis  with  the  large  mir¬ 
ror  is  given  by 

L'  =  (z^+x'2)4  (A-9) 


where 


z'  =  D'+2d(l-cos0) 

x'  =  2d  sin0 

Again,  for  small  0  ,  Eq.(A-9)  can  be  written 

L'  =  (D'2+4d202)^  ( A-10) 

Using  Eq.(A-8)  and  Eq.(A-7) 

L 1  =  [ (D-3d62 ) 2+4d2  0  2  ] ^  (A-ll) 

This  becomes 

L’  =  [D2+4d202+3d62(3d62-2D)p  (A-12) 

For  small  6  ,  using  the  binomial  approximation,  Eq.(A-12) 

becomes 

L'=  L+iAS( AS-2D)  (A-13) 
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Then 


AL  -  $AS( AS-2D)  (A-14) 

This  can  be  neglected  in  phase  terms  if 

ALk<<|  ( A-15 ) 

where  k  is  the  wave  number. 

Recalling  that  D  equals  d  in  the  case  being  considered, 
and  neglecting  terms  of  order  S'*  ,  it  is  found  that  the 

condition 

J 

;>,/  "  \  ( A-16) 

is  required  so  that  significant  phase  changes  do  not  occur. 

For  optical  frequencies  this  limits  6  to  very  small  values. 

Parenthetically,  a  tilted  off-axis  resonator  is  equiv¬ 
alent  to  an  off-center  non-confocal  resonator  (Ref  8:2242). 

A  tilt  can  then  be  viewed  as  a  rotation  about  the  center  of 
curvature  and  a  shift  of  the  feedback  mirror  along  the  op¬ 
tical  axis. 
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Appendix  B 


Virtual  Images  of  Feedback  Mirror  Edges 

The  field  amplitude  functions  f(z,x)  and  g(p,0)  in 
the  laser  cavity  are  expanded  in  terms  of  geometrical  ampli¬ 
tude  functions  and  edge  diffraction  amplitude  functions.  The 
diffraction  amplitude  functions  can  be  viewed  as  cylindrical 
wavelets  originating  at  the  virtual  images  of  the  feedback 
mirror  edges.  The  location  of  the  virtual  image  sources  for 
the  cylindrical  wavelets  can  be  determined  by  matrix  optics 
techniques . 

Given  the  geometry  of  the  off-center  unstable  confocal 
resonator  in  Fig.  (7),  the  image  of  any  point  (d,x)  on  the 
feedback  mirror  can  be  found  by  the  matrix  equation 

Je  T2RTj  (B-l ) 

where 

x  «*  transverse  coordinate  of  object  point 

6  15  angle  of  ray  leaving  object  point  with  respect  to 
optical  axis 

Tj  »  translation  matrix  from  object  point  to  mirror 

Rm  ■  reflection  matrix  for  imaging  mirror 

T2  ■  translation  matrix  from  mirror  for  exit  ray 
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Fig.  7  RAY  FROM  FEEDBACK  MIRROR 

x'  =  transverse  coordinate  for  exit  ray 
0'  =  angle  of  exit  ray  with  respect  to  optical  axis 
Substituting  for  the  translation  and  reflection  matrices 
(Ref  5:84-100)  yields 

x' 

e1 

where 

L  ■  translation  distance  parallel  to  optical  axis  from 
object  point  to  mirror 

L*  *  translation  distance  from  mirror  along  optical  axis 
R  ■  radius  of  curvature  of  mirror 


1  L' 
0  1 


1  0 
-f  1 


1  L 
0  1 


(B-2) 


75 


The  sign  convention  for  the  radius  of  curvature  is  chosen  to 
be  positive  for  concave  mirrors  and  negative  for  convex  mir¬ 
rors  . 


Performing  the  matrix  multiplication  of  T 
and  T2  in  Eq.(B-2)  yields 


X* 

1- 

A  R 

l  2l;l  +l'1 

X 

6’ 

2 

“  R 

i-  2k 

R 

e 

. . 

R  , 


(B-3) 


or 


“ 

- 
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X1 

Mu 

Ml  2 

X 

0’ 

M2  i 

M2  2 

e 

- 

- 

- 

.  . 

(B-4 ) 


The  matrix  for  an  image  forming  system  has  the  Mj2 
element  equal  to  zero. 

This  implies 

L  -  +  L'  =  0  (B-5) 


Given  L  and  R  ,  L'  can  be  determined.  The  distance  of 
L'  is  the  position  along  the  optical  axis  of  the  image  point 
from  the  reflecting  mirror.  For  L'  negative  the  image  is 
virtual  and  located  behind  the  reflecting  surface,  for  L* 
positive  the  image  is  real  and  in  front  of  the  mirror. 

For  an  image  forming  system  the  matrix  element  in 

Eq.(B-4)  is  the  transverse  magnification  MT: 
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- up 


MT  =  1  -  — -  (B-6 ) 

The  transverse  distance  of  the  image  from  the  optical  axis 
can  be  determined  from 


x*  =  Mtx  (B-7 ) 

So  the  object  point  (L,x)  will  be  imaged  at  (L',x')  rel¬ 
ative  to  the  imaging  mirror. 

For  the  off-center  unstable  confocal  resonator  geometry 
of  Fig. (7),  the  feedback  mirror  has  a  radius  of  curvature 

Rx  =  -2d  (B-8) 

The  large  mirror  M  has  a  radius  of  curvature 

R2  =  2(D+d)  (B-9) 

The  separation  between  the  mirrors  is  distance  D  .  The  feed¬ 
back  mirror  is  distance  d  from  the  focal  point. 

To  locate  points  relative  to  the  focal  point,  the  par¬ 
axial  coordinates  (p,x)  are  used  where  positive  p  is  the 
distance  to  the  left  of  the  focal  point  and  positive  x  is  the 
transverse  distance  above  the  optical  axis.  In  this  coordinate 
system,  a  point  on  the  feedback  mirror  is  (d,x) 

Given  an  object  point  (p0bj»xobj^  translation  dis¬ 

tance  Li  from  the  feedback  mirror  is 

L,  -  cobj-d  (B-10) 
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Then  using  Eq.(B-lO)  in  Eq.(B-5)  yields 


Lid 

L7+cT 


(B-ll) 


The  p  coordinate  of  the  image  point  is  then 


pim 


d+L' 


(B-12 ) 


The  transverse  magnification  of  the  feedback  mirror  using 
Eq. (B-8)  in  Eq.(B-6)  is 


MT 


d+L 


JL  s 


'im 

d 


(B-13 ) 


Using  this  in  Eq.(B-7)  yields 


d+L 


im 


obj 


pim 

d  xobj 


(B-14) 


Thus,  the  image  location  formed  by  the  feedback  mirror  is 
given  by 


*pim»xim* 


,  ,.T  •  d+L 
(d+Li,  — -  x 


obj 


(B-15) 


Similarly,  given  the  object  point  ^P0bj,xobj^  to  be 
imaged  by  the  large  mirror  M2 

L2  “  D+d-pobj  "  dM"Pobj  (B“16> 


.  »  _  L,  CD+d)  _  L,dM 
"  Lj -(D+d )  =  L2 -dM 


(B-17) 


78 


P,  *  D+d-L2  =  dM-L2 


(B-18) 


M  -  D+d-L2  pim 
T2  *  D+d  "  dM 


(B-19) 


D+d-L2 

D+d 


IAS  x 

dM  obj 


(B-20) 


where 


(B-21) 


The  object  points  on  the  feedback  mirror  undergo  transverse 
magnification  M  given  above. 

Then,  the  location  of  images  formed  by  the  large  mirror 
is  given  by 


*  <D+d-L-  W1  *obj> 


(B-22) 


‘  (dM-L>'  T1) 


(B-23) 


Using  Eq. (B- 16) through  Eq.(B-20),  the  virtual  image  of 
a  point  (d,x)  on  the  feedback  mirror  is  found  to  be 


(PiiXi)  -  (dM2 ,  xM) 


(B-24 ) 


The  image  of  this  point  (Pi,Xi)  produced  by  the  feed¬ 
back  mirror,  using  Eq.(B-lO)  through  Eq.(B-15)  is  located  at 
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(B-25 ) 


Cp_x  )  *=  (dM“*,xM'1) 

It  can  easily  be  shown  that  continued  reflection  of  the 
succeeding  virtual  images  yields  the  general  form 

( pn > xn )  =  (dM2n,xMn)  ;  n=±l,±2, . . . ,±N  (B-26) 

where  the  positive  subscript  indicates  the  nth  reflection  by 
the  large  mirror,  and  the  negative  subscript  is  the  nth  re¬ 
flection  by  the  feedback  mirror. 

Given  the  edges  of  the  feedback  mirror  located  at  (d,b) 
and  (d,-a),  the  virtual  images  of  the  edges  are  given  by 


(W  = 

(dM2D,bMn)  ;  n=±l,±2 ,  .  .  ,±N 

(B-27) 

<p0  »bo )  = 

(<J,b) 

(B-28) 

and 

(pn»“an>  *  (dM*n,-aMn)  ;  n=± 1 , ±2 , . . . ±N  (B-29) 

(Po»-ao)  "  (d»-a>  (B-30) 

It  can  be  seen  that  these  points  lie  on  the  halves  of  two  dif¬ 
ferent  parabolas  centered  on  the  optical  axis  as  shown  in  'Pig. 
(8).  These  parabolas  have  a  common  point  at  the  focal  point 
of  the  resonator.  This  is  the  (pn>bn)  and  (pn,-an)  point 
for  n -*■-<» 
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Appendix  C 


Asymptotic  Expansion  of  Integrals 
by  Method  of  Stationary  Phase 

This  appendix  outlines  the  development  of  the  asympto¬ 
tic  expansion  about  the  stationary  phase  point  of  integrals 
with  rapidly  oscillating  exponentials.  This  development  gen¬ 
erally  follows  that  of  Erdeli  (Ref  1:51-56).  The  general  form 
of  the  expansion  is  then  applied  to  the  integral  form  whose 
specific  asymptotic  expansion  is  used  in  Chapter  III. 


General  Expansion 

The  general  form  of  the  integral  to  be  expanded  is 


I(x) 


g(t)e 


ixh(t ) 


dt 


(C-l) 


where  x  is  a  large  real  variable  and  h(t)  is  a  real  func¬ 
tion  of  a  real  variable.  The  stationary  phase  point  is  defin¬ 
ed  as  the  point  t=t0  where  h'(t)=0  is  satisfied.  The  maj¬ 
or  contribution  to  the  integral  occurs  about  this  point.  For 
other  values  of  t  the  exponential  is  rapidly  oscillating 
and  the  average  contribution  is  very  small.  The  other  major 
contributors  to  the  integral  are  the  endpoints.  The  integral 
of  Eq.(C-l)  is  broken  up  into  two  integrals  at  the  stationary 
phase  point  t0  .  This  is  done  to  place  the  stationary  point 
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at  an  endpoint  of  the  integration,  this  is  needed  for  the  ex¬ 
pansion.  Then  Eq.(C-l)  becomes 

rt  0  r  8" 

I(x)  =  /  +  /  g(t)elxh(t)dt  (C-2 ) 

■  ot  to. 

Each  of  these  integrals  must  meet  four  requirements  be¬ 
fore  being  expanded.  First,  t0  is  the  only  stationary 
point  on  the  integration  interval.  Second,  h(t)  is  strict¬ 
ly  increasing  on  the  integration  interval.  Third,  the  end¬ 
points  of  the  integral  are  stationary  point  of  order  P-1 
for  the  lower  endpoint  and  a-l  for  the  upper  endpoint  where 
p  and  a  are  greater  than  or  equal  to  one.  These  orders 
can  be  found  by  expressing  the  derivative  of  h(t)  in  the  form 

h'(t)  =  ( t-a) p~ 1 ( 6-t )°~ 1H(t )  (C-3 ) 

where  H(t)  is  any  function  that  permits  h'(t)  to  be  written 
in  this  form.  The  fourth  requirement  is  that  g(t)  not  be 
zero  at  t=t0 

From  the  third  requirement,  integrals  of  the  form  of  Eq. 
(C-l)  fall  into  three  catagories: 

1.  For  a  stationary  point  at  the  upper  limit  of  inte¬ 
gration,  P=1  and  a>l 

2.  For  a  stationary  point  at  the  lower  limit  of  inte¬ 
gration,  p>l  and  a=l 
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3.  For  no  stationary  points  on  the  interval  of  inte¬ 
gration,  p=l  and  o=l 

An  integral  of  the  Eq.(C-l)  form  can  be  expressed 


g(t)eixh(t)dt 


B(x)-A(x) 


(C-4) 


where 


B(x)  =  /  [ 1-vC t ) ]  g(t)eixh(t)dt 

a+n 


(C-5) 


-A(x) 


Q-n 


v(t)g(t)eixh(t)dt 


(C-6) 


and 


v(t) 


1  ;a<t£a+ri 
0  ;  B-ri^t^B 


(C-7) 


The  v(t)  function  is  a  Van  der  Corput  neutralizer,  an 
N-times  differentiable  function  on  the  interval  of  integra¬ 
tion.  It  is  used  to  isolate  the  contribution  to  the  integral 
from  an  endpoint.  This  neutralizer  is  illustrated  in  Fig. (9). 

The  B(x)  term  is  evaluated  using  integration  by  parts. 

A  change  of  variable  is  made  where  the  new  variable  of  inte¬ 
gration,  w  ,  is  given  by 

h(t)  =  h(B)-a>°  (C-8 ) 


84 


From  this  the  new  limits  of  integration  become,  for  t=S 


o)j  =  0 


(C-9) 


and  for  t=(a+n) 


CO  2 


h(6)-h  (a+n)p 


(C-10) 


A  new  neutralizer  can  be  written 

6  ( to )  =  l-v(t) 


Its  values  at  the  endpoints  are 


,  U)=OJ] 


(C-ll) 


(C-12  ) 
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Also  the  function  y (o» )  is  defined 


y<“>  -  hs 


(C-13) 


Using  Eq.(C-8)  through  Eq.(C-13)  the  B(x)  integral  becomes 


B(x)  =  -e 


ikh( 6) 


/<*>2 

y(“ 


)  6(a))e_ixa)  da) 


(C-14 ) 


To  integrate  by  parts  the  following  are  chosen 


u  =  y(u)6(w)  ;  du  =  [y(a))6(u))]  (C-15) 


dv  =  e~ixilJ  ;  v  =  E^-x.io) 


(C-16) 


where  E  (-x,a))  is  some  function  satisfying  Eq.(C-16). 
The  integration  by  parts  of  Eq.(C-14)  yields 


B(x)  =  eixh(6)y(o))5(a))E  (-x,  a))eiXu>° 


0)2 


+e 


ixh(g) 


r  o)2 

)l  Ea(-x,a,)  ~  [y ( 0)) 6( 0))  ]  da> 
J  0 


(C-17) 


Here  it  should  be  noted  that  to  carry  the  asymptotic  expan¬ 
sion  to  higher  order  terms  the  integral  on  the  right  side  of 
Eq.(C-17)  is  integrated  by  parts.  For  this  analysis,  terms 
from  this  second  integration  are  neglected. 

Using  Eq.(C-12)  the  neutralizer  isolated  the  contribution 
to  the  B(x)  to  the  coi^O  endpoint  (from  t=8): 
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B(x)  -  eixh(6)y(0)E  f-x,0) 


(C-18) 


A  function  E0(-x,u>)  is  needed  which  satisfied 


h  *<,<-*•“>  -  8'ix“‘ 


(C-19) 


A  function  satisfying  this  is  given  by  the  contour  integral 


,  , Nn+i  /  ^  .  o 

)| — /  (z-w)neixz  dz 


(C-20) 


where 


y  -  exp  |  i  arg(z-u)  “  exP  |  -  55  | 


(C-21 ) 


Then  for  the  n®0  case 


0  -ixz 


(C-22 ) 


d  .  ,  v  -ixwc 

♦-!<“>  "  “e 


(C-23) 


The  E  ( — x . u> )  function  is  then 
0  ’ 


E^C-x.u)  - 


(C-24 ) 


To  get  an  expression  for  Eo(-x,0)  in  Eq.(C-18)  the  con¬ 
tour  integration  of  6_n_ t C )  must  be  performed  along  the 
path  from  u>  to  infinity  at  angle  yu  to  the  real  axis 
which  is  illustrated  in  Fig. (10). 

A  change  of  variable  is  made,  letting 
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Fig.  10  PATHS  OF  CONTOUR  INTEGRATION 


z  -  te 


—  1 TT 

w~ 


With  this  change  of  variable  and  letting  w-0 
becomes 

*  (0)  -  e 

♦_n_jVO)  ■  -e  I  t  e  dt 

Jo 

Letting  s“-xt°  ,  this  equation  is  then 


♦_n_,(°)  ■  -ox 


f 


X  1  s 

’« 


e~sds 


Eq, 


(C-25) 

(C-20) 

(C-26) 

(C-27) 
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The  Integral  term  in  this  equation  is  a  gamma  function  as 
defined  by 


r(m) 


e_sds 


(C-28 ) 


where  m  -  .  Then  for  n  ■  0  and  using  Eq.(C-24) 

o 

Eo(-x,0)  ■  o”lx  °  e  2or(o) 

From  this,  for  a“l 


(C-29) 


iir 

E,(-x,0)  -  x“‘e  r(l)  ■ 


(C-30 ) 


And  for  o*2 


E2(-x,0) 
Similarly,  the 


-  -ix-V^  r<*) 


(C-31 ) 


-A(x)  term  of  Eq.(C-6)  can  be  found  to  be 


-A(x)  -  -eixh(a)k(0)Ep(x,0)  (C-32 ) 

where  Eq.(C-6)  has  been  integrated  by  parts  using 

h(t )  -  h(a)  +  up  (C-33) 

k(u)  -  g(t)  ^  (C-34) 

and  also  using  the  following  to  get  the  final  form 
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(C-35  ) 


Ep(x,u)  -  ♦_l (u) 

Ep(x,0)  -  -p_1x  P  e  ^  r(p/  (C-36 ) 

Then  for  p-1,2 

E»(x'0)  "  IS  (C-37) 


E,(x,0)  -  -J  nr 


(C-38 ) 


Finally,  Eq.(C-4)  yields  the  general  form  of  the  asympto¬ 
tic  expansion  about  the  stationary  point,  t0  ,  using  Eq. 
(C-18)  and  Eq.(C-32): 


re 

I  g(t)eixh(t)dt  -  y(0)Eo(-x,0)eixh(8) 


(C-39) 


-k(O)Ep(x,O)eixh(0l) 


An  Application  of  the  Asymptotic  Expansion 

Given  an  integral  of  the  form  of  Eq.(C-l)  that  is  to  be 
approximated  by  an  asymptotic  expansion  about  its  stationary 
phase  point,  five  steps  are  used: 

1.  The  stationary  point  of  h(t)  is  found  such  that 
h’(t0)-O 

2.  The  integral  is  broken  up  into  two  integrals  at  its 
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stationary  point,  t#  ,  as  Eq.(C-2). 

3.  Any  changes  of  variable  are  made  that  are  necessary 
to  ensure  h(t)  in  each  of  the  two  integrals  is 
strictly  increasing  over  the  interval  of  integration. 

4.  The  value  of  p  and  o  for  each  integral  is  deter¬ 
mined  by  expressing  h'(t)  in  the  form  of  Eq.(C-3). 

5.  Each  of  the  integrals  is  solved  as  in  Eq.(C-39). 

For  this  example,  the  integral  to  be  solved  is 


f6 

(x)  -  I  g(t)  exp  j  -ix(at *-2bt+c) |  dt 


. (C-40) 


The  stationary  point  is  determined  from 


h’(t)  -  2(at-b) 


( C-4 1 ) 


This  yields  the  stationary  point,  t „•  ^  ,  for  h'(t)-0 

It  is  assumed  that  a< 

Then  Eq.(C-40)  is  broken  up  into  two  integrals  about  the 
stationary  point  — 


g(t)  exp 


-ix(at*-2bt+c) J  dt 

(C-42 ) 


For  the  second  integral  h(t)  is  strictly  decreasing  on  the 
interval  of  integration.  A  change  of  variables,  t--t  ,  is 
performed  to  make  this  strictly  increasing.  Then  Eq.(42)  can 
be  expressed  as  I«I|+It  ,  where 
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(C-43) 


r/a  i  , 

■  I  g(t)  exp  J  -ix(atl-2bt+c) |  dt 

rb/.  .  . 

[2  ■/  g(-t)  exp  j  -ix(ati+2bt+c)  |  dt  (C-44) 


The  first  of  the  two  integrals,  ,  is  now  expanded. 

In  the  form  of  Eq.(C-l),  h(t)  is  identified  as 


h(t )  -  -(at*-2bt+c)  (C-45) 

Expressing  h'(t)  in  the  form  given  in  Eq.(C-3)  yields 

h'(t)  -  2a(£  -t)  (C-46) 

From  this  it  can  be  seen  that  p«l  ,  am 2  ,  and  H(t)»2a 

Then  using  Eq.(C-39) 


Ii  -  y(0)E2(-x,0)eixh'a/ 
-k(0)E,(x,0)eixh(a) 


(C-47) 


From  Eq. (C-45) 

-(!)  -  r-< 


(C-48) 


h(a)  ■  -(aa*-2ba+c) 


(C-49) 


Next,  the  y(0)  term  is  determined.  From  Eq.(C-8)  it  is 
found  that 
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(C-50) 


bi 


a 


at 


■2bt  + 


b* 

a 


Then,  using  the  quadratic  formula 


(C-51) 


The  root  is  chosen  so  that  ^  causes  the  resultant  y(0) 
term  to  be  positive.  In  this  case  the  positive  root  is 
chosen  and  so 


dt  „  _1 


(C-52  ) 


Using  Eq.(C-13),  Eq.(C-51)  and  this  expression,  the  y(<*>) 
term  at  w«0  is 

y(0)  -  g(|)  (C-53) 

/r 


Next  the  k(0)  term  is  evaluated.  From  Eq.(C-33)  (Re¬ 
calling  that  p-1  for  this  case),  Eq.(C-45),  and  Eq.(C-49) 

u  -  -at *  +2bt+( aa* -2ba )  (C-54) 


Using  the  quadratic  formula 


t  -  |?(a,-2|a+|i-H)i  (C-55) 

Here  the  root  is  chosen  so  that  ^  will  cause  k(0)  to  be 
negative.  For  this  case  the  positive  root  is  chosen,  then 
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"  -J(a2a2-2baa+b2-au)”^ 


(C-56) 


Using  Eq. (C-34) , Eq. (C-55) ,  and  this  expression,  the  k(u) 


for  u-0  becomes 


k(0)  -  ig(a )  aa-b }”* 1 


(C-57) 


Now  Ij  from  Eq.(C-47)  can  be  expressed,  using  Eq.(C-31), 
Eq. (C-37) ,  Eq. (C-48) ,  Eq.(C-49),  Eq.(C-53),  and  Eq.(C-57) 


*>  ■  5g(!)  e*p|  -ix(c-li)l 

+g(o) |2ix(aa-b) ]  exp  j  -ix(aa2-2ba+c)J 


(C-58) 


By  a  similar  procedure  I2  is  found  to  be 


*'(!)  ViF  exp  |-ixH!)l 

+g(0)[-2ix(a6-b)] 

X  exp  J  -ix(ae*-2be+c) J 


(C-59) 


Finally,  the  integral  from  Eq.(C-40)  is  expressed  in  an 
asymptotic  expansion  to  terms  of  order  x”1  as 


f  8 

I  g(t)  exp  j -ix(at2-2bt+c)|  dt 

■  (ir)4  *(!) exp  I  -lx(c-l-)l 
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+g(a)[2ix(aa-b)]  exp  |  -ix(aa2-2ba+c)| 
+g(8)[-2ix(a$-b)J  exp  | -ix(aB2-2be+C) J 
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Appendix  D 


Asymptotic  Approximation  of  Diffraction  Integra.!, 
for  Off-Center  Unstable  Resonator 

This  appendix  develops  the  asymptotic  expansion  for  the 
diffraction  integral  given  in  Eq.(39)  where  f(D,an)  is  re¬ 
placed  by  the  form  in  Eq.(26).  The  form  of  Eq.(39)  is  then 

g(p,0)  =  Cl  =  C(Ii+I2+I3)  (D-l) 


where 


ik 

"3T 


( an-bn ) 2 


X  rn(d,an)  exp 


I 

2 


t  <V*>“ 


n»l 


exp 


ik  la^aj,2 
^  “a-* 


X  sQ(d,an)  exp 


(D-2  ) 


(D-3) 


(D-4  ) 
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+F(nj)  [ 2 ± ( Arj 2  — B ) J  exp  J  i(An22~2Bn2+C)J 


Note  that  this  agrees  with  the  asymptotic  expansion  form 
Eq.(C-60)  developed  in  the  previous  appendix  by  another  meth¬ 
od. 

Approximation  for  I, 

When  I,  is  put  into  the  form  Eq.(D-7)  the  following 
terms  are  found 


n,  *  -1  ,  n2  =  - 


(D-9) 


F(n)  =  rn(d,an) 


(D-10) 


Ai  "  2irFea  -d.(P-d)+P^pn-d) 
(p-d)(pn~d) 


(D-ll) 


B,  -  2ttF 
i  ea 


p0d(pn-d)+dbn(p-d) 


a(p-d)(pQ-d) 


(D-12 ) 


c.  -  w-«) 


Using  Eq.(D-ll)  through  Eq.(D-13)  it  can  be  shown  that 


B L 
A, 


d  bn(p-d)+p0(pn-d) 
a  d(p-d)+p(pn-d) 


(D-14) 
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i 


7» 


l(d,an)  oxp 


-ik 

T~ 


hi 

(FI) 


(D-5 ) 


An  asymptotic  series  approximation  of  integrals  , 

I2  ,  and  Ij  can  be  made  when  they  are  in  the  form  (Ref  3: 
1541) 

I  -  /*n*  F(n)  exp  | -it(A'n2-2B’n+C' )|  dn  (D-6) 

*n» 

whore  t  is  large  valued.  In  Ij  ,  I2  ,  and  Ij  let 
t  correspond  to  -2nF  whoro  F  is  largo  valued  and 
given  by  Eq.(3). 

The  Eq . (D-6 )  can  be  written 


I  -  rnz  F(n)  exp|i(An,-2Bn+C)|  dn 
•'ni 


(D-7) 


whoro 


A  -  -2ttF  A1 
ea 


B 

C 


-2ttF  B' 
ea 

-2ttF  C  ' 
ea 


The  asymptotic  expansion  of  this  integral  when  terms  of  high¬ 
er  order  than  F  ”l  are  neglected  is 

ea 


[nrl  F(f)  oxp  I  i(c"  f“)l  (D-8) 

+F(rn )  |-2i( An  ,-B)  ]  'exp  |  i(An1a-2Bn1+C)| 
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t 


/a 


f(d,an)  exp 


-ik  (e-  d  n)' 

2  tnr 


(D-5) 


An  asymptotic  series  approximation  of  integrals  Ij  , 
I2  ,  and  I 3  can  be  made  when  they  are  in  the  form  (Ref  3: 
1541) 

I  =  /*nz  F(n)  exp  [  -it(A'n2-2B'n+C'  )|  dn  (D-6) 
n  j 

where  t  is  large  valued.  In  lx  ,  I2  ,  and  I3  let 
t  correspond  to  -2nF  where  F  is  large  valued  and 
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given  by  Eq.(3). 

The  Eq.(D-6)  can  be  written 


I 


F(n)  exp  i(An2-2Bn+C) 


dn 


(D-7) 


where 

A  =  -2rrF  A' 
ea 

B  =  -2t:F  B' 

ea 

C  *  -2ttF  e  • 
ea 

The  asymptotic  expansion  of  this  integral  when  terms  of  high¬ 
er  order  than  F  “l  are  neglected  is 
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ItI  F(!)  eHi(c-!1)l  (D-8) 

+F(m)  [-2i(Anj-B)]  'exp  (  i(Ani2-2Bni+C)| 
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! 


+F(n2 )  [2i(An2-B)]  lexp  |  i(An22-2Bn2+C)| 

Note  that  this  agrees  with  the  asymptotic  expansion  form 
Eq.(C-60)  developed  in  the  previous  appendix  by  another  meth 
od. 

Approximation  for  I, 

When  Ij  is  put  into  the  form  Eq.(D-7)  the  following 
terms  are  found 

Hi  -  -1  n2  -  |  (D-9) 


F(n)  =  rn(d,an) 


(D-10) 


2itF  d(p-d)+p(pw-d) 

6a  <P-d)(PQ-d) 


(D-ll) 


B,  =  2irF 
i  ea 


p0d(p  -d)+db_(p-d) 


n 


a(p-d)(pn-d) 


(D-12) 


C.  -  2"Fea  <D'13> 


Using  Eq.(D-ll)  through  Eq.(D-13)  it  can  be  shown  that 


El 

Ai 


d  bn(p-d)+p0(pn-d) 
a  d(p-d)+p(pn-d) 


(D-14) 


u 

Approximation  for  1 2 

The  limits  of  integration  are  the  same  as  Eq.(D-9).  When 
I2  is  put  into  the  form  of  Eq.(D-7)  the  following  terms  are 
found 


F(n)  ■  sn(d,an) 


*«  -  ™.a 

».  -  2"rea  p»d^rd? 


C2  -  2ttF  p9  d  (Pn-d)-t-danz(p-d) 
ea  a2(p-d)(pn-d) 


It  then  follows  that 


-  d  P9(pn-d)-an(p-d) 
*2  a  d(p-d)  +p(pn-d) 


c,‘  xh  ■  j  dfp-d)-dp<p -d> 


AjHi-Ba  “  - 


+  p(a+9d) 

Pn-d  p-d 


Aam-Ba  ■  jy— 


(D-21) 

(D-22 ) 

(D-23) 

(D-24 ) 

(D-25 ) 

(D-26) 

(D-27) 

(D-28) 
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( 


Aan»l 

-2B2Hi+Ca  “  j 

(a_an 

)2  +  fi(a+0d) 

k-3 

T  d(T-*) 

Aa  n  2  * 

-2Ban2+C2  “  j 

(>>♦»„ 

)2  .  p(b-0d) 

l*v* 

r  d(p-d) 

(D-29) 


(D-30) 


Substituting  into  Eq,(D-8)  yields 


I* 


i^-'W 

n-l  *  l 


2tt i  d(p-d)(pn-d) 


,  * 


kaJ  d(p-d)+p(pn-d) 


X  sn  (“•« 

x  expl2h  I 


♦  sn(d,-a) 


*  S»j 
x  «p  I  £  (imtL  ♦  I 


-i 


+  sn(d,b) 


Mf^S  +  SW) 


-- 1 


(D-31) 


Approximation  for  I» 

The  limits  of  integration  for  I»  are  given  in  Eq.(D-9). 


When  I)  is  put  into  the  form  of  Eq.(D-7)  the  following  terms 
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are  found 


F(n)  -  f(d,an) 

(D-32) 

A*  “  2ltFea 

(D-33) 

B>  ’  2"Fea 

(D-34) 

C,  .  M  e»L£- 

ea  a2 (p-d) 

(D-35) 

Then  it  follows  that 

r  0d 

(D-36) 

C-|ji  -0 

(D-37) 

-  -5s  tfgtf2 

(D-38) 

A,n2-B,  -  5-  dl’p-'jy^ 

(D-39) 

A.ru'-aB.n.+C,  -  |  P4<p-aj 

(D-40) 

A.n.’-aB.n.+c,  -  $ 

(D-41) 

Substituting  into  Eq.(D-8)  yields 

r,.4  -  h 


When  Eq. (D-2) ,  Eq.(D-20),  Eq.(D-31),  and  Eq.(D-42)  are  sub¬ 
stituted  into  Eq.(D-3)  the  expression  for  g(p,8)  near  the 
feedback  mirror  is  obtained.  The  final  result  is 


x  - 1  Is  *  $*&■)  i 


♦  r„<d’b>  lika(^9  + 

+  Yj  (pn"d)“*(llj£  dipldivp^p^l 

n»l  \ 

*  °*H  r  ,C&---v  I 


sn(d,-a) 

IMS^s 

.  p(a+0d)\  1 
d(p-d)  /  | 

“P|?( 

la~an.ll  + 

pn-d 

fiW)l 

sn(d,b)  [ 

ik“(^9 

♦»r 

e„„  1  ik  / 

(b+a„)2  . 

pfb-0d)2\ I  \ 

exp  1  5"  ( 

pn-d 

S(p-B  )l  J 

+  ^jla5:  p  (p-d)  1  *(d’6d) 

♦  *«*.-»>  i-s^r 
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Appendix  E 


Elements  of  Coefficient  Matrix  A(u) 


This  appendix  illustrates  several  of  the  elements  of  the 
coefficient  matrix  A(y)  of  Eq.(140).  The  coefficients  are 
of  a  very  complex  form  but  their  values  could  be  calculated  on 
a  computer  since  the  eigenvalues,  \i  ,  are  the  only  unknown 
parameters. 

From  Eq.(l^l),  the  A(jj)  matrix  is  square  and  of  dimen¬ 
sion  (2N+1)  .  It  can  be  expressed  as 


A  - 


1  .  X 


1  »  1 


1  •  2 


2  »  2 


1  ,N+1 


*M1,2N  M1,2N+1 


•  •  • 

M2N+1,1  M2N+1 ,  £ '  *  ^2N+1 ,  N+i 


2N+1.2N+1 


(E-l) 


where  U.  _  is  a  matrix  element. 
k,n 

Similarly,  the  q  vector  can  be  expressed 


q(xN) 


q(x, ) 
q(x0) 

qt*.,) 


q(x-N+l) 

«(x-N> 


(S-a> 
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From  Eq. (139)  several  examples  of  coefficient  matrix 
elements  are 


v,  -  CE<yww>'N-1 

(E-3) 

"lit  CEN-1AN-1(x+N)uN 

(E-4) 

N+l 

M1,N+1  2CE°1-p 

(E-5) 

Ml,N+2  "  CE-1A-1(x+N)u 

(E-6) 

“l,2N+l  "  CE-NA-N(x+N)uN 

(E-7) 

M2,2“  CEN-1AN-1(xN-1)ijN~1-1 

(E-8) 

M2N+1,1  "  CENAN(x-N)yN 

(E-9) 

'2N+1.2N+1  "  CE-NA-N(x-N)uN_1 

(E-10) 

Then,  using  Eq.(134)  through  Eq.(138),  several  examples  of 
these  matrix  elements  in  more  detailed  form  are 
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